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Abstract 

Massive stars live fast and die young. They shine furiously for a few million 
years, during which time they synthesize most of the heavy elements in the 
universe in their cores. They end by blowing themselves up in a powerful 
explosion known as a supernova (SN). During this process, the core collapses 
to a neutron star or a black hole, while the outer layers are expelled with 
velocities of thousands of kilometers per second. The resulting fireworks often 
outshine the entire host galaxy for many weeks. 

The explosion energy is eventually radiated away, but powering of the new- 
bom nebula continues by radioactive isotopes synthesized in the explosion. 
The ejecta are now quite transparent, and we can see the material produced 
in the deep interiors of the star. To interpret the observations, detailed spec- 
tral modeling is needed. This thesis aims to develop and apply state-of-the-art 
computational tools for interpreting and modeling SN observations in the neb- 
ular phase. This requires calculation of the physical conditions throughout the 
nebula, including non-thermal processes from the radioactivity, thermal and 
statistical equilibrium, as well as radiative transport. The inclusion of multi- 
line radiative transfer, which we compute with a Monte Carlo technique, rep- 
resents one of the major advancements presented in this thesis. 

On February 23 1987, the first SN observable by the naked eye since 1604 
exploded, SN 1987A. Its proximity has allowed unprecedented observations, 
which in turn have lead to significant advancements in our understanding 
of SN explosions. As a first application of our model, we analyze the ^^Ti- 
powered phase (t > 5 years) of SN 1987A. We find that a magnetic field is 
present in the nebula, trapping the positrons that provide the energy input, and 
resulting in strong iron lines in the spectrum. We determine the ^^Ti mass to 
1.5+1] J X 10"^ Mq. From the near-infrared spectrum at an age of 19 years, 
we identify strong emission lines from explosively synthesized metals such 
as silicon, calcium, and iron. We use integral-field spectroscopy to construct 
three-dimensional maps of the ejecta, showing a morphology suggesting an 
asymmetric explosion. 

The model is then applied to the close-by and well-observed Type IIP SN 
2004et, analyzing its ultraviolet to mid-infrared evolution. Based on its Mg 
I] 4571 A, Na I 5890, 5896 A, [O I] 6300, 6364 A, and [Ne II] 12.81 urn 
nebular emission lines, we determine its progenitor mass to be around 15 Mq. 
We confirm that silicate dust, SiO, and CO have formed in the ejecta. 

Finally, the major optical emission lines in a sample of Type IIP SNe are 
analyzed. We find that most spectral regions in Type HP SNe are dominated by 
emission from the massive hydrogen envelope, which explains the relatively 
small variation seen in the sample. We also show that the similar line profiles 
seen from all elements suggest extensive mixing occurring in most hydrogen- 
rich SNe. 
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1 Introduction 




Supernovae (SNe) are explosions of stars that have ceased their fusion pro- 
cesses. They represent astrophysical laboratories where a large and diverse 
set of physical processes take place. By studying them, we can learn about 
the formation of neutron stars and black holes, test theories of stellar evolu- 
tion and nucleosynthesis, and understand the chemical enrichment histories of 
galaxies. We can determine their contribution to the dust and cosmic rays in 
the universe, and refine their use as standard candles for cosmological distance 
measurements. 

SNe emit radiation at all wavelengths, from radio to gamma-rays. The di- 
verse set of observations is matched by an equally diverse set of models and 
theories needed for their interpretation. Initially, the expanding fireball is hot 
and opaque, and models for this phase use radiation hydrodynamics to simu- 
late the coupled evolution of radiation and matter, providing information on 
the explosion energy, the ejecta mass, and the radius of the progenitor. Scat- 
tering models may be applied to the outer atmospheric layers to determine 
their density profiles and composition. Later on, as the SN expands and cools, 
its inner parts become visible, glowing from input by radioactive elements 
produced in the explosion. The matter is now further from thermodynamic 
equilibrium, and models have to consider a large set of physical processes to 
compute the state of the gas. In this phase, we can determine the morphology 
of the SN ejecta, the amount of elements produced by the nucleosynthesis, the 
amount of mixing and fragmentation that has occured in the explosion, and 
possibly observe signatures of the compact object formed at the center. 
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Introduction 



This thesis deals with modehng of such late phases. The research field is 
still in its infancy, as SNe are dim and hard to observe at late times. For many 
years, the early phases of the nearby SN 1987A have been the only observa- 
tions studied in depth. But due to recent advancements in observational sur- 
veys, the number of SNe with spectral data in the later phases is increasing. 
Combined with developments in atomic data calculations and easy-access par- 
allel computing, the field is moving rapidly forward. We are now coming to a 
point where accurate spectra from realistic explosion models can be computed 
and compared with an increasing number of well-observed events, allowing 
constraints to be placed on which type of stars explode as which type of SNe. 
In the coming years, the development of several robotic wide-field surveys 
will further increase the inflow of data, and theorists will have their hands full 
to analyze them. 

The specific aims of this thesis have been to develop state-of-the-art com- 
putational models for nebular-phase SNe, and apply these models to derive 
properties of hydrogen-rich core-collapse SNe. We initially investigate the 
very late phases of the famous SN 1987A, which has mainly been analysed 
in its early phases before. We then investigate a recently obtained sample of 
other Type IIP spectra, with special focus on the interesting object SN 2004et. 

To elucidate the nature of SNe, we have to understand several subjects in 
physics and astronomy. Chapter 2 deals with the evolution of the massive stars 
that are their progenitors, providing the underpinning for what SN ejecta ac- 
tually contain. While the study of SN may eventually come to revise stellar 
evolution theory, these are the current ideas and the perceived uncertainties 
in them. Chapter 3 provides a broad overview of SNe, including a histori- 
cal perspective and observational aspects such as classification. Chapters 4 
and 5 contain a rather detailed description of the fundamental theory incorpo- 
rated into the computational model developed. In Chapter 4, 1 describe how 
the physical conditions in the SN are determined by a balance between ra- 
dioactive, thermal, and radiative processes. Chapter 5 describes the treatment 
of the radiation field, which represents a major advancement in the modeling 
compared to earlier models developed in Stockholm and elsewhere. Finally, 
Chapter 6 provides a summary of the results that we have obtained by apph- 
cation of the model. 

Most things in the universe evolve slowly. SNe are different, providing 
entertainment on a time-scale more suitable for the human mind. 
They evolve into something new in about the same time it takes 
you to write a paper about what they were before. I hope to have cap- 
tured some of their fascination in this thesis, and that you will enjoy reading it. 

A.R.J. 

Stockholm, November I, 2011 
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2 Evolution of massive stars 



'Astronomy? Impossible to understand and madness to investigate!' 

Sophocles, 420 B.C. 

Stars are gravitationally confined fusion reactors. Inside them, light elements 
are converted into heavier ones, releasing energy that provides pressure sup- 
port against the gravity and makes the star shine. Low-mass stars like the sun 
take a few billion years to deplete their supply of hydrogen, after which grav- 
itational contraction occurs until helium fusion begins. During this burning 
process, helium is converted to oxygen and carbon. In the next contraction 
stage, the densities become so high that the star stabilizes due to electron de- 
generacy pressure, and a white dwarf is formed. 

Massive stars (M > 8 M©), on the other hand, never reach that stabilizing 
point. They ignite also their central supplies of carbon and oxygen, and 
after that also heavier elements until the core has been transformed to iron. 
The final structure of the star is onion-like, with layers made up by the 
ashes of the various burning stages (Fig. 12.1b . In this chapter, I review 
the major aspects of the lives of massive stars u p until the forma t ion of 



the ir on core. Specific st e llar e volution models (jHashimoto et al 



Hirschi et al. (2004, H04). iNomoto et all (.1997. N97), Woosle v et al 



H89), [W ooslev & Weavei] (Il995[ WW9 5), 'Thielemann et al.' (Il996 



(11989 . 
T96) 

tool . 



WHW02V lWooslev & Hegeil (120071 WHO?')^ are occasionally referenced. 



2.1 Fusion 

Fusion reactions are mediated by the nuclear force, which has an effective 
range of /? ~ 10^'^ cm. Nuclei must therefore come within this distance to 
each other for reactions to occur. Such close encounters are opposed by the 
repulsive electric force between the positively charged nuclei. We may attempt 
to estimate the critical temperature required for fusion to occur by equating the 
kinetic energy to the electrostatic potential energy, and solving for T: 

Ekin ~ Epot ^keT^ , (2. 1) 

K 

where Z is the proton number, q is the charge unit, and kg is Boltzmann's 
constant. For hydrogen (Z = 1), we obtain T ~ 10^" K. At the beginning of 
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Figure 2.1: Schematic final structure of a massive star that has synthesized elements 
up to iron (Fe). 

the 20th century, calculations showed that stars were not this hot inside (the 
sun is only ~ 10^ K at its center), and fusion was therefore rejected as their 
power source. However, two effects make fusion effective even at ~ 10^ times 
lower temperature; some nuclei have much higher kinetic energies than kT , 
and secondly quantum mechanical tunneling allows reactions to occur also 
for sub-energetic particles. It was therefore not until the 1930s, when quan- 
tum mechanics was thoroughly understood, that fusion was established as the 
power source of stars. 



2.1.1 Hydrogen burning 

The net effect of hydrogen fusion, starting at T ~ 10^ K, is 

4^H^'^He. (2.2) 

In addition, two positrons and two neutrinos are created. The path taken is 
not four protons simultaneously colUding into one helium nucleus, but rather 
a sequence of successive buildups. In massive stars, the dominant path is the 
CNO cycle, in which primordial abundances of carbon, nitrogen, and oxygen 
are used as catalysts. The main branch is 

'2C(p,7)i3N(e+v,)"C(p,7)i4N(p,7)i5o(e+v,)i5N(p,4He)i2c . (2.3) 

For a complete cycle, ^^C is recycled in the last step. However, while the 
cycle is running, the abundances of the involved nuclei are altered. The rea- 
son is that the various reactions occur at different rates, with the ^^N(p, yY^O 
reaction being significantly slower than the other ones. In equilibrium, the 
number of ''^N nuclei is therefore enhanced and the number of ^^C nuclei is 
suppressed. Another reaction occurring is ^^0(/?, y) ^^F(e+ v^) ^^0(/?,'* He) '^N, 



2.1 Fusion 
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which converts O to N. The net result is that the carbon and oxygen abun- 
dances are suppressed during CNO burning, whereas the nitrogen abundance 
is enhanced. These shifts in abundances are retained when fusion reactions 
eventually cease. 

Table 12.1] shows the composition of the H-burning ashes from a few differ- 
ent models. We can see the alteration of the original C, N, and O abundances. 
This alteration explains why oxygen can be a coolant in the outer (unbumt) 
hydrogen envelopes of SNe, but not in the He zones (paper IV). 



El. 


WH07-15 


WHW02-15 


H04-20 


He 


0.98 


0.98 


0.98 


C 


1.9 X 10-4 (0.08) 


n/a 


2.5 X 10-4 (0.10) 


N 


9.0 X 10"^ (11) 


0.011 (14) 


n/a 





1.8 X 10-4 (0.03) 


n/a 


3.0 X 10-4 (0.05) 


H 


-10-^ 


n/a 


n/a 



Table 2.1: Composition (mass fractions) of the hydrogen-burning ashes from a few 
stellar evolution models. The mass of the star is given after the reference. The change 
relative to solar abundances for the CNO elements are given in parenthesis. All other 
elements have their natal mass fractions. 



During the hydrogen-burning phase, the star is said to be on the main se- 
quence. Massive stars will belong to classes O or B, with surface temperatures 
of 20,000-40,000 K, and luminosities of 10^ - 10^ Lq. The time-scale for 
this phase is about 



tMS ~ 40 ( ) million years . (2.4) 



2.1.2 Helium burning 

As hydrogen is eventually exhausted in the core, the star contracts and heats 
up. At r ~ 10^ K, helium fusion begins, occurring by the triple-a reaction, 

3 4He ^ ^^C , (2.5) 

and also by 

i2c + 4He^i60, (2.6) 

which explains why '^C and ^^O are the major carbon and oxygen isotopes 
on earth and in the universe. The rate of the last reaction is still only known 
within a factor two or so, and is one of the major uncertainties in stellar evo- 
lution models, both for the carbon/oxygen yields, and for the yields of heavier 
nuclei produced from them in later burning stages (WHW02, WH07). Table 
12.21 shows the composition of the C/0 zone formed by helium burning from a 
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few different models. It is from this zone we can see emission from molecular 
CO in SNe (paper III). 



El. 


H89-20 


WH07-12 


H04-20 





0.76 


0.79 


0.65 


C 


0.21 


0.18 


0.30 


Ne 


0.015 


0.011 


5.0 X 10-^ 



Table 2.2: Composition (mass fractions) of the helium-burning ashes, taken from afrw 
stellar evolution models ( only the three most common elements are listed). The mass 
of the star is given after the reference. 



During the helium-burning phase, the stellar envelope expands and cools, 
and the star becomes a red supergiant (RSG), with photospheric temperature 
on the order of a few thousand degrees. 

2.1.3 Carbon burning 

At r > 8 X 10^ K, carbon nuclei fuse as 

2 ^ ^^Mg (2.7) 

^ 20Ne + a (2.8) 

^ 23i^^_^p (2.9) 

^^Mg + n. (2.10) 



The neutrons, protons, and a-particles will cause further reactions, resulting in 
a large set of nuclides. The main products are neon and magnesium, comple- 
mented by some carbon, sodium, silicon, and aluminium. The zone is referred 
to as the 0/Ne/Mg zone later in the thesis, and we shall see in paper III how 
this zone is key for linking SNe to their progenitors. 

Table l2.3l shows the composition of this zone from a few different models. A 
major down-revision in modern calculations is the amount of silicon produced 
in this zone (compare the H89 value with the other ones). Since silicon is an 
efficient coolant (see e.g Kozma & Fransson 199 8a,), this has major impact on 
the thermal emission from this zone. One may also note that this is the main 
production site o f neon and sod ium in the universe, so one may compare the 
solar ratio of 33 ("Loddersl booi) to the model outputs. 

From the carbon burning stage and on, the burning time-scales shorten sig- 
nificantly compared to the hydrogen and helium-burning time-scales. The rea- 
son is that internal energy is now efficiently converted to neutrinos that escape 
freely. 
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El 


H89-20 


WH07-12 


H04-20 


N97-18 


o 


0.73 


0.66 


0.58 


0.66 


Ne 


0.12 


0.25 


0.31 


0.21 


Ms 

o 


0.075 


0.043 


n/a 


0.037 


Na 


7.0 X 10""^ 


7.7 X 10-3 


n/a 


0.010 


C 


5.7 X 10-3 


3.9 X 10-3 


4.5 X 10-3 


0.045 


Al 


8.0 X 10-3 


4.0 X 10-3 


n/a 


5.0 X 10-3 


Si 


0.020 


3.0 X 10-3 


4.0 X 10-3 


2.5 X 10-3 



Table 2.3; Composition (mass fractions) of the carbon-burning ashes, from a few stel- 
lar evolution models ( only the three most common elements are listed). The mass of 
the star is given after the reference. 



2.1.4 Neon burning 

Before the oxygen ignition temperature is reached, neon nuclei start photo- 
disintegrating to oxygen, as the photon energies are now in the MeV range. 
The released a-particles are captured by other neon nuclei to make magne- 
sium, silicon, and sulphur. The reactions begin at T ~ 1.7 x 10*^ K, and can be 
summarized as 

20Ne-F7 ^ + a (2.11) 
20Ne-Fa 24]y[g (2.12) 
24Mg-fa ^ 28si (-2.13) 
2^Si-Fa ^ 32s _ (2.14) 

(2.15) 

The net effect is that two neon nuclei are converted to one oxygen nucleus and 
one magnesium nucleus, with some of the magnesium nuclei being further 
converted to silicon and sulphur. The composition after this burning stage is 
mainly oxygen, silicon, sulphur, and some magnesium, and is referred to as 
the 0/Si/S zone in the thesis. Table 12.41 shows the composition in this zone 
from a few models. 



2.1.5 Oxygen burning 

At r ~ 2. 1 X 10^ K, oxygen bums as 

2i6o ^ 32s (2.16) 



31s + n (2.17) 

3ip-Fp (2.18) 

28Si + a (2.19) 

3<'P + d. (2.20) 
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El. 


N97-18 


H04-20 


WH07-15 





0.80 


0.73 


0.82 


Si 


0.10 


0.08 


0.10 


S 


0.025 


0.01 


0.015 


Mg 


0.032 


n/a 


0.044 



Table 2.4: Composition (mass fractions) of the neon-burning ashes, from a few stellar 
evolution models ( only the three most common elements are listed). The mass of the 
star is given after the reference. 



The main products are silicon and sulphur (an Si/S zone), with traces of chlo- 
rine, argon, potassium, calcium and phosphorus. From this stage on, the burn- 
ing ashes will be further modified in the eventual SN explosion, so I do not list 
the chemical compositions. The nucleosynthesis abundances of these ashes 
are strongly non-solar, and the fact that they are not observed supports the 
previous staternent th at these ashes will be further modified before ejection 
dWooslev etal.lll972h . 



2.1.6 Silicon and sulphur burning 

After oxygen burning is complete, all further reactions occur by photodis- 
integrations followed by a-captures (similar to the neon burning). Above 
T ~ 3.7 X 10^ K, silicon and sulphur nuclei fuse with a-particles to produce 
iron-group nuclei through steps of ^^Ar, ^°Ca, '*'*Ti, ^^Cr, ^^Fe and ^^Ni. The 
burning occurs under conditions close to nuclear statistical equilibrium (NSE), 
where all electromagnetic and strong reactions are in detailed balance. The 
main isotope produced under NSE is the one that is most tightly bound for the 
given neutron exces^. Typically, the neutron excess is close to zero, and that 
nucleus is then ^^Ni. 



2.1.7 Neutron and proton capture 

During the star's evolution, nucleosynthesis occurs partly via the slow (s) pro- 
cess, in which nuclei capture free neutrons. By a series of such captures, el- 
ements more massive than ^^Ni can also be made. The neutrons are mainly 
produced by successive a-captures on ^'^N nuclei, and the s-process occurs 
mainly at the end of helium burning (WHW02). 

Protons are less likely to fuse with nuclei due to the electric repulsion. How- 
ever, some important reactions still occur, one example being proton capture 
by ^^Ne to produce ^^Na. This isotope is radioactive with T = 3.75 years, and 



The neutron excess is the number of neutrons minus the number of protons, divided by the sum 
of them. 
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serves as a source of radioactivity in the 0/Ne/Mg zone (paper II). Another 
example is proton capture by ^^Mg to produce ^^AI. 

An important consequence of the neutron and proton captures is that the 
abundances of elements not participating in the main fusion cycles themselves 
may be altered from their natal values. For example, in the 0/Ne/Mg zone the 
nickel and cobalt abundances are strongly enhanced, with significant amounts 
of the radioactive isotope ^"Co produced (paper II). 



2.2 Mass loss 



Observations of massive stars indicate t hat they are los ing material from their 
surfaces (e.g.. IWright & Barlow! 1 19751 lLamerslll98lh . Their optical spectra 
often show P-Cygni profiles (blueshifted absorption and redshifted emission), 
which arise only from material in rapidly outflowing motion. From the line 
profiles, gas velocities higher than the escape velocity of the star can be in- 
ferred. Another method to detect the outflowing gas is by its radio free-free 
emission, or by far-infrared (FIR) emission from dust forming in the outflow. 

By studying a large sample of stars, the dependency of the mass-loss 
rate on stellar parameters such as luminosity, temperature, and metallicity 
can be determined. Such relations have been derived both for main 
sequence stars and for RS G s. A commonly used relation is the one by 
Nieuwenhuijzen & de Jager (Il99(3l) . A high dust production efficiency, 
combined with low surface gravities, makes mass loss rates higher for RSGs 
than for main sequence stars. 

Theoretically, mass loss on the main sequence is understood as 
radiation pressure acting on atomic lines. This idea was first proposed by 
Lucv & SolomonI (179 70 ). and subsequently quantitative models have been 
developed that accurately reproduce the observed mass loss rates. The 
importance of metallicity is obvious for this type of mass loss. 

In the RSG phase, understanding of the mass loss mechanism is poorer. 
Here, radiation pressure probably acts on dust grains forming in the cool enve- 
lope. It is also possible that large-scale convective motions and pulsation s play 
a role, as recent imaging of Betelgeuse suggests (IKervella et al.ll201 ih . Self- 
consistent modeling is challenging, and existing models still have to make a 
large number of assumptions. 

More massive stars lose more mass, in total, than less massive ones 
(INieuwenhuijzen & de Jageii Il990h . Models suggest that the final envelope 
mass is therefore quite insensitive to the ZAMS mass of the star, being 
6 — 10 Mq for non-rotating stars below 25 Mq, and steadily smaller for 
more massive ones (H04, WH07). Rotation increases the mass loss rate, and 
always leads to a smaller final mass (H04). 

For stars more massive than ~35 M0, also the helium core may be affected 
by mass loss (WH07). A 100 M0 star ends up with a final mass of only ~6 
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M0, with no hydrogen or heUum. In the models of WH07, the heUum core 
reaches its maximum size (13 Mq) for a 45 Mq star. We therefore find the 
important result that most types of (single) massive stars reach final masses of 
5 — 15 Mq at the end of their evolution. 



2.3 Convection 

If stars were stable objects, it would be relatively easy to compute their evo- 
lution through the different burning stages. Unfortunately, they are usually 
not. Convection refers to the turbulent mass motions that represent the sec- 
ond way for the star to transport energy besides radiation. A detailed theory is 
so far not available, because the solutions to the equations of hydrodynamics 
become chaotic for turbulent flows, preventing meaningful results to be ob- 
tained. However, phenomenological models such as the mixing-length theory 
are available which are believed give reasonable results. 

Convection transports energy and alters the chemical composition by mix- 
ing gas between different layers. H-burning ashes can be convectively dredged 
up into the hydrogen envelope, altering the fractions of helium and CNO el- 
ements (Sect. [2?TTT] ). Similarly, He-burning ashes can be dredged up into the 
He zone, altering its content of oxygen, carbon, and neon (Sect. 12.1.2b . This 
mixing typically does not reach all the way through the He layer, but produces 
a He/C zone in its lower parts. In paper IV, we see a strong [C I] 8727 A line 
in the spectrum of SN 2008bk, that appears to arise from such a dredged-up 
layer. 

Large uncertainties currently exist in the treatment of convection. Two dif- 
ferent stability criteria for convection may be derived, the Schwarzs child crite- 
rion and the Ledoux criterion. These are identical for a perfect gas of constant 
composition. But if radiation pressure is important, or if there is a compo- 
sition gradient, the Schwarzschild criterion may predict instability whereas 
the Ledoux criterion may not. In this situation, the region is said to be i^em/- 
convective. It is no t known exactly how the gas behaves in this situation. 



Woosley & Weaveri (119881) allow diffusion for the composition, but not for 
the thermal energy, a recipe followed in many models also today. 

Another source of uncertainty regards the so called convective 
over-shooting. The stability criteria define exact radii in the star where the 
convection-driving buoyancy force disappears. But when the convective 
bubbles reach these points, they have non-zero velocities and will continue 
a bit passed the limit. The treatment of t his over-shooting c ritically 



influences all aspects of the evolution of the star ( SchaUer et al Jl 1 992 ) . Strong 



over-shooti ng leads to a more massive c ore, a faster burning rate, and a higher 



luminosity (IMowlavi & Forestinilll994l) . Unfortunately, the predicted extent 



of over-shooting ranges from being neglegible to almost doubling the size of 
the convective regions. 
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2.4 Rotation 

It is well known that mas sive stars on the main sequence rotate rapidly 



(WHW02). iFukudal (Il982h find typical rotation velocities of 200 kms"^ 
Rotation affects the star mainly during its hydrogen and helium burning 
phases, with the later burning stages transpiring too quickly for the rotational 
instabilities to have any impact (H04). Rotation leads to larger helium cores 
for a given stellar mass and metallicity, caused by an enhanced mixing of fuel 
into the nuclear burning core. For example, a 15 M0 star has a helium core of 
3.0 Mq in the non-rota ting models of WHW02, but of 4.9 Mq in the rotating 



models of iHeger et al. Idooo). For a 20 Mq star, the numbers are 5.0 M0 and 



7.8 M0. Rotation can also induce outward mixing of burning ashes to the 
surface, and increase the mass loss rate both by the centrifugal effect and by 
increasing the core luminosity (WHW02). 

All models used in this thesis are non-rotating ones. When we connect the 
amount of metals observed in the SN to the progenitor mass of the star (Paper 
III), we therefore actually derive an upper limit. We do know, however, that 
rotation does not have any strong influence on the chemical composition of 
the nuclear burning ashes (H04), so the metal line ratios will not change much 
by computing spectra from rotating progenitors. 

We may finally note that rotation influences the final point in the HR dia- 
gram. H04 finds the SN progenitor of a 20 M0 star to be red, yellow, and blue, 
for rotation velocities of 100, 200, and 300 kms^', respectively. 



2.5 The final structure of the star 

An important property of the sequence of fusion reactions described in Sect. 
12.11 is that each ash is only partly burned in the next burning stage. The star 
is therefore develops an onion-like structure, with successively heavier ele- 
ments further in. The details of the final structure depends on the specific 
prescription used for rotation, convection, mass loss, and uncertain reaction 
rates. Currently, the uncertainties associated with convection (Sect. 1231 ) is the 
largest sources of diversity in stellar evolution models (WHW02). 

Fig. 12.21 shows the final structure of a 20 Mq non-rotating star as computed 
by H04. Note the stratification arising from the nuclear burning sequence H— >■ 
He, He^C/0, C^O/Ne/Mg, Ne^ 0/Si/S, Si/S, and Si^ ^^Ni, modified 
by convective dredge-up of carbon into the He zone, and helium into the H 
envelope. 



2.6 Concluding remarks 

Today, we know that all atomic nuclei in the universe, apart from hydrogen, 
helium, lithium, beryllium, and boron, are created inside stars. Up until 1957, 
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Figure 2.2: Final structure (mass fractions) of a 20 M© non-rotating star From 



Hirschi et al. (2004). 



this was only one of four theories being considered for the origin of the el- 
ements. By the key role played by helium nuclei in nucleosynthesis (Sect. 
12.1b . we have already come a long way to understand the large abundances of 
elements with even proton numbers in the universe. 

In Chapter 2, 1 describe how the elements produced are ejected into the in- 
terstellar medium, to eventually become building blocks for new stars, planets, 
and humans. 
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'In the second year of the Chung-phing reign period, the 10th month, on a Keui-day 
(December 7, AD 185), a guest star appeared in the midst of the Southern Gate. It 
seemed to be as large as half a bamboo mat, and showed five colours in turn, now 
brightening, now dimming. It diminished in brightness little by little and finally 
disappeared in the sixth month of the Hou-year (Hou-nian, 24 July to 23 

August AD 187): 

Astrological Annals of the Hou Han Shu ( China) 

Before the invention of the telescope at the beginning of the 17th century, not 
much seemed to happen in the night sky. Only two types of objects provided 
the occasional drama; comets and guest stars. 

Historical records of guest stars date almost two millenia back in time. The 
oldest one is from year 185, when Chinese astronomers recorded a bright new 
star that visited the night sky for a few weeks (the excerpt above). Similar 
events happened in years 386 and 393. The next records are from Chinese and 
Arabic astronomers in 1006 and 1054. These new stars were taken to convey 
important messages, surely altering the fates of many individuals and maybe 
even civilizations. 

Today, we know that these guest stars are novae and supernovae (SN). The 
Swedish astronomer Knut Lundmark established, in 1920, that SNe are about 
a thousand times intrinsically brighter than the novae. Later it has been clar- 
ified that novae are surface eruptions on stars, while SNe are the complete 
explosive disruptions of them. 

That SNe are explosions of stars was first suggested by Baade and Zwicky 
in 1939, who asked themselves what would happen to a massive star once it 
reached the degenerate iron core stage described in the previous chapter. In- 
spired by the discovery of the neutron by James Chadwick in 1932, their an- 
swer was that the core would collapse to a neutron star, with the gravitational 
binding energy released in that process ejecting the outer layers in a violent 
explosion. This type of event is called a core-collapse supernova. There is 
also another class, the thermonuclear supernovae, which are the explosions of 
carbon-oxygen white dwarfs. 

Today, we use our telescopes to zoom in on the locations of the historical 
guest stars mentioned above. Indeed, Baade and Zwicky were right; we find 
the expanding debris of violent explosions, sometimes with a pulsating neu- 
tron star at the center! Fig 13. II shows what is seen in the direction of the year 
185 and 1054 guest stars. 
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Figure 3.1: The remnants of SN 185 and SN 1054 (the Crab nebula). Image credit : 
Chandra/NASA. 



3.1 Collapse and explosion 

If fusion reactions were to cease inside the Sun, it would slowly start contract- 
ing as thermal energy diffused out, lowering the pressure support inside. In the 
core of a massive star, thermal energy is not lost by such slow radiative dif- 
fusion, but by immediate neutrino losses. Once fusion stops, things therefore 
happen quickly. Apart from the neutrino losses, which follow pair annihila- 
tions and electron captures, internal energy is also lost by endoergic photo- 
disintegrations of the iron nuclei. In less than one second, ten million years 
of fusion reactions are undone in the core as the iron nuclei are disintegrated 
into a-particles. In its own immense gravitational field, the core collapses on 
a free-fall time-scale of ~ 100 milliseconds. Matter crashes towards the center 
at about a quarter of the speed of light, compressing the core from the size of 
the earth to the size of a city. At that point, nuclear densities are reached, and 
the repulsive part of the nuclear force kicks in. The inner parts of the collaps- 
ing core come to a screeching halt, forming a proto-neutron star (PNS). The 
sudden deceleration launches a shock wave 20 km from the center that travels 
upstream through the infalling remainder of the iron core, which accretes onto 
the PNS after being shocked. The shock loses energy as it fights its way out 
through the inflow, but copious neutrino emission from the PNS pushes it on. 
For a few seconds, the PNS emits as much energy in neutrinos as the rest of 
the universe does in light! The shock wave runs through the rest of the star, 
exploding it. Approximately 10^^ ergs of energy is endowed to the outer lay- 
ers, which are expelled at thousands of kilometers per second. A supernova is 
born. 

The details of the explosion mechanism outlined above are still being 
worked out. No spherically symmetric models succeed in achieving the 
required energy transfer from the neutrinos to the shock wave, instead 
showing the whole star collapsing to a black hole, simply disappearing from 
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the night sky. It has therefore become clear that some sort of multidimensional 
effects must be crucial. Several types have been proposed, including 
rotationally-driven explosions along the star's rotation axis, convection, and 
acoustic oscillations. In the last few years, 2D and 3D models have emerged 
that show that so called Standing Accre tion Shock Instabilities (SASI) occur 



in the first few seconds after explosion (iBIondin & Shawll2007l IScheck et al. 



20081). These pick out one or a few directions by spontaneous symmetry 



breaking, enhancing the neutrino deposition efficiency by convective motions 
along these modes. In Paper I, the 3D mapping of the ejecta of SN 1987A 
shows a morphology consistent with SASI-like instability. 



3.1.1 Explosive nucleosynthesis 

The shock wave running through the star is strongly supersonic and radiation- 
dominated. One may show that such a shock accelerates material to 6/7 of the 
shock velocity, and heats it to (Sedoy 1959) 

r=4300pi/V//^ K, (3.1) 
where Vg is the shock velocity, and p is the densi ty, which is of ord er lO' g 



cm^^ in the layers outside the collapsed core (e.g. iHirschi et al.ll2004t) . Com- 
bining this value with the typical shock velocity ~ 10^ kms~\ we can 
estimate T ~ 8 x 10^ K. Since this is higher than the the silicon-fusion limit 
(Tsi/s ~ 4 X 10^ K), the shock burns the material to iron-group elements. The 
burning occurs close to NSE (Sect. [2?L6l) . producing ^^Ni as the main isotope. 
However, the density soon falls to a regime where so called a-rich freeze-out 
produces large amounts of '^He nuclei, as well as radioactive isotopes such 
as ^^Ti. The excess of a-particles is built up as as the inefficient triple-alpha 
reaction falls out of equilibrium. In paper II, we discuss the signatures of this 
explosive burning in the late-time spectrum of SN 1987A. 

It can be shown that the shock decelerates if the density falls off less steeply 
than a r^^ power law, but accelerates otherwise. As the shock travels through 
the star, both decelerations and accelerations occur. However, for all reason- 
able density profiles, the quantity p(r)'/'^Vv(r)'/^ in Eq. l3.1l steadilv decreases, 
so the post-shock temperature decreases. Explosive burning occurs in reverse 
order compared to the hydrostatic burning. For a 20 M0 star, the whole sil- 
icon/sulphur zone, and part of the oxygen zone, is burnt under NSE to ^^Ni, 
followed by partial burning of the oxygen zone to S i/S (oxygen burning) 



and 0/Mg/Si/S (neon burning) ( Hashimoto et al. 19891) . The carbon burning 



regime occurs in regions where only trace carbon exists, and by the time the 
shock reaches the helium and hydrogen layers, the shock temperature is too 
low to burn them. In total, elements heavier than magnesium are mostly made 
explosively, and lighter ones mostly by hydrostatic burning. 

Since the shock is radiation-dominated, the energy density is given by e = 
ar^, and we can estimate how much material will be burned by approximating 
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the temperature as constant and equating 



(3.2) 



where Es is the shock energy, a is the radiation constant, and R is the radius. 
For Eg = 10^' ergs and T = T^i/s, the solution is R = 3700 km, which roughly 
corresponds to the radiu s within which all material will be burned to iron- 
group elements (Thielem ann et al.lll996h . For typical pre-SN models, the cor- 
responding mass coordinate is 1.7 M0. With the mass cu^ for the neutron 
star being around 1.6 Mq, the estimate for the amount of ^^Ni ejected is then 
0.1 M0, in good agreement with the observed values. Similarly, the explosive 
oxygen-burning and neon-burning occurs out to ~6400 km and ~ 12000 km, 
with corresponding mas s coordinates ~1.8 and ~2.1 Mq, for a 20 Mq star 
(IThielemann et al.lll996h . The empirical correlation between explosion energy 
and ^^Ni mass, shown in paper IV, can likely be understood from these sim- 
ple considerations. If the mass cut stays roughly constant, a higher explosion 
energy simply keeps the shock hot enough for silicon-burning over a larger 
volume. 



3.1.2 Reverse shock and mixing 

Several types of hydrodynamical instabilities may occur in the expanding gas. 
The Rayleigh-Taylor (RT) instability occurs when lighter parcels of gas are 
accelerated into denser ones. The interface between the parcels then fragment 
and the fluids get mixed, resulting in finger-like structures. The instability 
arises w hen the pre s sure gradient an d the density g radien t have opposite di- 
rections. IC hevalieij (fl97 6) and C hevalier & KleinI (Il978h showed that such 



conditions indeed arise in the flow of hydrogen-rich SNe. iHerant & Woosley 
d 19941) showed that the RT mixing is stronger in RSGs compared to blue su- 
pergiants. The mixing in SNe with RSG progenitors (paper III, IV), should 
therefore be at least as strong as observed for SN 1987A, which had a blue 
progenitor. 

These early simulations used one-dimensional codes until the shock 
had reached the hydrogen envelope, and only then mapped the ejecta to 
a multi-dimensional grid, necessitating that perturbations be put in by 
hand. Recently, however, computational advancements have allowed the 
multi-dimensional computations to cover also the early phases. As mentioned 
in Sect. 13.11 these simulatio ns find strong inst abilitie s arising from the 



earliest times. Kifonidis et al. 1(|2003) and Kifonidis et al .l (I2OO6) compute the 



evolution in 2D, finding the SASI instability to arise already during the first 
second, leading to a strongly asymmetric shock. These asymmetries provide 
seeds for the RT instabilities at the Si/O and 0/He interfaces, fragmenting 



'The division line between material falling back to form the compact object, and material that 
is ejected in the SN. 
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and mixing them within the first few minutes after core bounce. The He/H 
interface is, however, found to be mixed mainly by Richtmyer-Meshkov 
instabilities, which arise when a shock front impacts a composition interface 
at an oblique angle. These results are important in showing that also the core 
zone interfaces are subject to instabilities, which leads us to apply strong 
artificial mixing throughout the core region of the input models used in this 
thesis. 

One may show that the diffusion time scale for mixing on the microscopic 
(i.e., atomic) scale , is long compared to the hydrodynamic timescale 
(IFryxell et al.lll99ll) . In addition, models assuming that microscopic mixing 
does not occur generally agree better with observations than models tha t 
assume that it does (e.g. Fransson & Chevalier 19891 Maeda et al. 2007 ). 
Any significant microscopic mixing between the metal zones and the helium 
and hydrogen zones can be ruled out based on the empirically established 
production of CO and SiO, as helium ions destroy the se molecules 
dLe pp et all [l99nl iLiu & Dalgamol Il996l lGearhartetaDll999b . The mixing 
discussed above is therefore believed to occur only on macroscopic scales. 
The difference is crucial for spectral modeling, since even small amounts of 
certain elements may completely change the spectrum formed in a particular 
zone. 

As the ejecta expand, the internal energy steadily decreases through ra- 
diative and adiabatic losse s. Internal energy is, h owever, also resupplied by 
radioactive decay of ^^Ni. iHerant & Bend (1199 ih find that the ^^Ni clumps 
expand on a time-scale of ~1 week, inflating themselves to hot, low-density 
bubbles. This result, and others, leads us to assume a large filling factor for 
the iron clumps in all models in this thesis. 

After some time, pressure forces become neglegible and hydrodynamical 
interaction ceases. Each parcel of gas then continues coasting on whatever 
trajectory it is on. The SN is then said to have entered the homologous phase. 



3.2 The diffusion phase 

As the shock breaks out through the surface of the star, the SN becomes visi- 
ble. The initial burst will be hot and peak in X-rays, but soon the gas expands 
and cools to emit mostly ultraviolet and visible light. Only one SN has been 



caugh t in the act of the actual shock breakout, SN 2008D (ISoderberg et al. 



20081). Th e early UV/optical coo ling phase has been observed in several ob- 



jects (e.g. lStritzi nger et al.l l2002l) . 

The highly ionized nebula is initially optically thick, and the energy in- 
side can only slowly diffuse towards the surface. The spectrum is that of a 
blackbody scattered by lines in the atmosphere. The major difference to the 
spectrum of an ordinary (hot) star is that the absorption and emission lines are 
much broader due to the high expansion velocities, which completely domi- 
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nate the thermal velocities. Many lines also form P-Cygni profiles, where the 
velocity field causes a blue-shifted absorption through combined with emis- 
sion on the red side. 

The early light-curve of the SN evolves in a way that depends on how much 
internal energy has been deposited by the shock passage, what fraction of 
this energy will be radiated rather than lost to adiabatic expansion, and on 
what time-scale that emission will occur. The major dete rminant for the latter 
two q uantities is the compactness of the star, M/R (e.g. 'imshennik & P opov 



1992r) . For a RSG, about 1% of the internal energy will be radiated away, or 



~ 10'*^ ergs. If the time-scale for emission is one month, the luminosity is 
then ~ 10^ L0, still matching that of an entire galaxy. The radioactive ele- 
ments produced may also influence the diffusion phase light-curve, especially 
for compact stars where the explosion energy is rapidly lost to adiabatic ex- 
pansion. 



3.3 The steady-state nebular phase 

After a few months, the expanding nebula becomes optically thin in the con- 
tinuurrH, and is said to have entered the nebular phase. The remaining internal 



Continuum absorption processes are those tiiat occur for a wide range of photon energies, such 
as electron scattering, photoionization, and free-free absorption. See Chapter 5 for more details. 
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energy can then be efficiently radiated away, and if no other energy source is 
present the SN will rapidly fade away to undetectability. Fortunately, the ra- 
dioactive elements produced in the explosion provide a new power source. 
As we saw in Sect. 13.1.11 the initial shock burning produces ^^Ni, which is 
radioactive with a life-time of 8.8 days. Its daughter nucleus, ^^Co, is also 
radioactive with a life-time of 111.5 days. Most SNe enter the nebular phase 
when the ^^Ni has decayed and ^^Co is the dominant power source. It provides 
a powering of, for t ^ 8.8 days: 

L(56Co) = 3.4 X lO^Le ( ^^^1 e-'l''' ' ^ , (3.3) 
y (J. 1 Ai0 J 

where M(^^Ni) is the mass of ^^Ni synthesized in the explosion, typically 
around 0.1 M0. The decay products thermalize to produce UVOIR 
(UV-optical-near-infrared) radiation from the SN. 

One can show that all atomic processes are fast in the early nebular phase 
(Chapter 4), which means that the SN almost instantaneously re-emits all the 
energy put in by radioactivity. We can therefore refer to this phase as the 
steady-state nebular phase. As long as the decay particles are trapped by the 
ejecta, the light-curve then follows the time-evolution of the ^^Co nuclide, 
L{t) oc d jj^ (j^jg phase, we do not need to know the history of the 

nebula to predict its output, just the instantaneous radioactive input. 

In general, the nebula will still be optically thick in various lines even in 
the nebular phase. This fact has been one of the key drivers for this thesis 
project; to develop models that include transfer through these lines throughout 
the ejecta. In paper II, we show that this line transfer indeed has a strong 
influence on the emerging spectrum. 

Can the line opacity trap radiation so that the output is actually lagging the 
radioactive input? The empirical fact that the bolometric light curves do not 
show any evidence of such time-delay suggests no. The velocity field makes 
transfer through individual lines a local process (see Sect. I5.2.21 ). so a photon 
in an optically thick line does not have to diffuse through the whole nebula. 
In Sect. 15.1.21 we show that lines are active over a length scale ~ 10^^^ of the 
size of the nebula, so the mean-free -path A is 

X ~ -^^^—11 , (3.4) 

where V is the expansion velocity scale, and T is the line optical depth. The 
average number of scatterings equals the optical depth T in the Sobolev ap- 
proximation (Sect. [5l2l2l ). so the time the photon is trapped in the line is 

fscatter = — = 10-'— <10-'? , (3.5) 

c c 

where we have used V /c ~ 10^^ in the last step. We see that the trapping time 
is independent of the optical depth of the line, and neglegible compared to 
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the evolutionary time. Line-to-line scattering will, however, introduce longer 
trapping times. A photon that scatters in N lines will, in worst case, be trapped 
for a time 



yt 



fscatter~A^— ~ IQ-^M , (3.6) 

c 

which can potentially violate the steady-state assumption if N is greater than 
a few. In practice, N is kept small by the fact that fluorescence occurs with 
a relatively high probability for most lines, t ransferring the photons to longer 



wavelengths where the line opacity is lower (IPinto & Eastman. 2000.) . 



3.4 The time-dependent nebular phase 

After several years, the radiative output from the SN no longer matches the 
instantaneous input by radioactivity, bec ause the reprocessing is going too 
slow ( Axelrod.l98a Iciavton et al.lll992l [pransson & Kozma 19931 The SN 
is then said to have entered the time-dependent nebular phase. From this time 
on, solutions of the time-dependent equations for temperature, ionization, and 
excitation are required, rather than their steady-state variants (see Chapter 4). 
This means that the history of the nebula now influences its output at any given 
time. 

Regions of low density and low ionization are the ones to be firstly affected 
by this time-dependence. The code developed in this thesis does not take time- 
dependence into account. In paper I and II, where the outer envelope of the 
SN has entered the time-dependent phase, we compute conditions in those 
layers with the time-dependent code by KF98. 



3.5 Supernova classes 

One can estimate that about one SN occurs per second somewhere in the uni- 
verse. Of these, we currently discover about one in 100,000, or a few hundred 
per year. The discovery rate is expected to increase dramatically in the coming 
years thanks to robotic wide-field surveys such as Pan-Starrs, The Supernova 
SkyMapper, Large Synoptic Survey Telescope (LSST), and the Palomar Tran- 
sient Factory (PTF). 

The empirical classification of SNe are divided into an initial branch of Type 
I (hydrogen lines present) and Type II (hydrogen lines absent). The Type I 
class then divides into Type la (strong Si II 615 nm line). Type lb (helium lines 
present), and Type Ic (helium lines absent). The Type II branch subdivides into 
Type IIP (a plateau in the light-curve). Type IIL (a linear decline of the light- 
curve), and Type Iln (narrow lines present). The Type la class corresponds 
to thermonuclear SNe, whereas all other classes correspond to core-collapse 
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Figure 3.3: Observed fractions of SN classes in a volume-l imited sample (left) and a 
magnitude-limited sample (right). Data from .Li et alJ (1201 ih . 



events. Some SNe not fitting into the standard classification scheme are given 
the peculiar (pec) designation. 

Fig. 13.31 shows the observed fractions of the various SN classes. Since 
brighter SNe are more easily detected, the observed fractions of those 
(magnitude-limited) is higher than the actual fraction (volume-limited). 



3.5.1 Type II 

The most common clas s in the univer se is the Type II class, which makes 



up ~57% of all events (iLi et all 120111) . Only in the last decade or two have 
decent nebular phase spectral data sets been obtained for Type lis, with SN 
1987A long being the only well-observed event. A listing of well-observed 
Type IIPs can be found in paper IV. Some of these appear to have similar 
^^Ni-masses and ejecta velocities as SN 1987A, but others seem to come from 
much weaker explosions, with low velocities and low ^^Ni-masses. 



3.5.2 Type Ib/c 

Explosions of stars that have lost their hydrogen envelope, or both the hydro- 
gen and helium envelopes, are called SN Type lb and Ic, respectively. Due to 
their lower fre quency of 19% (volume-limited) and 4% (magnitude-limited) 
dLi et al.11201 11), Type Ib/c SNe have even fewer nebular data sets than the IIPs. 
However, with their small envelopes and large metal cores, they represent an 
important class to analyse with respect to nucleosynthesis. 
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3.5.3 Type la 

Type la SNe arise as a white dwarf accretes matter until is reaches the Chan- 
drasekhar limit, at which time it ignites thermonuclear fusion of carbon and 
oxygen under degenerate conditions, leading to a run-away explosion. The las 
differ from core-collapse events in that they produce more ^^Ni (~ 1 M0 ver- 
sus ~0.1 M0), and reach higher ejecta velocities (~0.1c versus ~0.01c). 

This thesis is mainly focused on the core-collapse SN class. However, the 
code developed (described in Chapters 4 and 5, and in papers II and III) has 
been tested on a Type la explosion model and compared to the 'typical' Type 
la SN 2005cf. The spectrum is shown in Fig. 13.41 The general agreement is 
encouraging for the future application of the code for Type la modeling. 
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Figure 3.4: The spectrum of the Type la SN 2005cf at day 94 p ost-explosion (black 
line), compared to the model spectru m of explosion model W7 tNomoto et al.lll984t) 
(red line). From lMaurer et al.l (1201 lb . 



3.6 SN 1987A 

On February 23 1987, the first SN to be seen by the naked eye since 1604 
exploded in our satellite galaxy, the Large Magellanic Cloud. The proximity 
(50 kpc) gave an opportunity to study a stellar explosion in unprecedented 
detail, and SN 1987A continues to be carefully studied today. 

The early spectrum of the SN showed clear hydrogen lines, establishing it 
as a Type II event. However, the light-curve did not evolve as any other Type 
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II previously observed, being unusually dim and rapidly evolving. It became 
clear that the exploded star must have been more compact than a RSG, which 
are the usual progenitors of Type II SNe. 



3.6.1 The progenitor 

The progenitor star, S k —69° 202, had luck ily been registered in pre-explosion 
surveys of the region ( Gilmozzi et alJll987 ). the first time a SN progenitor had 
been caught on a photograph. It was confirmed to be a compact and hot star, 
a blue supergiant o f spectral class B3 la, with radius 50 Rq and luminos- 
ity 10^ Lq ( Arnett et all 1 1989) . It did not have any particularities associated 
with it (the region was actually observed regularly for over a century prior 
to 1987), confirming the theoretical expectations that the late burning stages 
would transpire too quickly to affect the appearance of the star. 

Stellar e volution codes w ere scrambled to see what models could explain Sk 
—69° 202. IWoosleyl (Il988h found a 20 Mq star with reduced semi-convection 
to end as a blue supergiant with the correct luminosity. Other models in- 
voked a binary system where a merger between a 15 — 20 Mq star and a 
~3 Mf7^ main-sequence star took place ^ 2 x 10^ yea rs prior to explosion 
dffillebrandt & Meyei][l989l IPodsiadlowski et al.lll990h . The exact nature of 
the progenitor system still remains unknown. For our modeling of SN 1987A 
in paper I and paper II, we assume that the progenitor had a mass of 19 M0. 
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3.6.2 Spectral analysis 

The analysis of SN 1987A had as of late 2011 spawned some 1100 refereed 
scientific articles, which comprises a body of information almost impossible 
to adequately overview. I mention here a few of the most important res ults. A 
good review of its first five years of evolution can be found in McCrayi (Il993h . 
Some overview of modeling results can also be found in paper II. 

Light curve 

Hydrodynamical modeling of the light-curve indicated a hydrogen envelope 
of size 7 — 10 MfT^ and an explosion energy of 0.8 — 1.5 x 10^^ e rgs (IWoosley 



1988L IShigevama & Nomotol ll990L iBethe & Pizzocherol Il990h . The shock- 
deposited energy dominated the light-curve for the first four we eks, but due 



to the rapid adiabatic cooling, radioactivity soon took over (e.g. lArnett et al 



19891) . From multi-band observations in the steady-state neb ular phase, the 
^^Ni mass was determined to 0.069 MgH dBouchet et al.lll99lh . 



Spectra 

Most emission lines in the nebular spectra showed similar lin e profiles, 
suggesting a similar velocity distribution (e.g. iMeikle et al.lll993h . Some of 
the clearly identifiable optical lines were Ha, [O I] AA6300, 6364, [Ca II] 
AA7291, 7323, [Ca II] IR-triplet, [Fe II] A7155, 7172, and Mg I] A4571. The 
absence of Hj8 in the spectrum could be explained by models showing that 
the optical depths f or the B aimer lines were high, resulting in conversion of 
Hj8 to Paa -i- Ha ( Xu et al.lll992 ). The hydrogen recombination line ratios 
in Type IIP SNe therefore follow the unusual 'Case C scenario, with thick 
Lyman and Balmer series. 

Not even the hydrogen lines showed any flat tops, indicating that hydro- 
gen must have been mixed down to almost zero velocities. At the same time, 
the iron lines showed broad wings suggesting that some of the ^^Ni had been 
mixed far out into the envelope. Re cently, the first 3D simulation of a realis- 
tic progenitor ( Hammer et al. 201(i) succeeded in reproducing the high-speed 
^^Ni bullets inferred to be present in the ejecta. 

Molecules 

For the first time, molecules were detected in a SN. Carbon monoxide (CO) 
was detected in its fundamental ba nd (4.6 jxva) and first o vertone band (2.3 
jUm) a few months after explosion (ISpyromilio et al.lll988h. The flux peaked 



at200 days and persisted to at least day 775 (IWooden et al.lll993l) . iLiu et al. 
(Il992h estimated a total CO mass of 1 x 10^^ Mq. This is in agreement with 



models showing that the non-thermal electron p opulation limits the condensa- 
tion efficiency to a few percent (.Clayton et al.li2001.) . This in turn means that 



^For D=50.1 kpc and color excess E(B-V) = 0.15 mag. 
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CO cannot be important in zones with low carbon abundance, and that the 
carbon in the C/0 zone is mostly available for dust formation. 

Silicon monoxide (SiO) w as detected in its fundamental band (8.1 /im) at 



day 160 (lAitken et al.lll988l) . persisting to day 519. Liu & Dalgarno (1996| ) 
estimate a mass of 5 x 10~^ M0, implying that neither SiO locks up much 
material. 

Molecules have since been detected also in other core-collapse SNe, carbon 
monoxide in 1995ad, 1998S, 1998dl, 1999em, 2002hh, 2004et, 2004fj, and 
2005af, and silicon monoxide in SN 2005af and 2004et. In paper III, we use 
these observations to motivate an assumption of CO and SiO cooling in the 
C/0 and 0/Si/S zones. 

Dust 

Several indicators showed that dust formed in the supernova during its second 
year, see paper II (Sect. 2.6) for a review. Dust formation has since been ob- 
served in several other SNe (1998S, 1999em, 2003gd, 2004dj, 2004et, 2006jc, 
and 2010jl), as well as in SN remnants Cas A, the Crab, Kepler and 1E0102.2- 
7219. 



3.6.3 Circumstellar interaction 

The presence of narrow emission lines in the spectra quickly made it clear that 
the star was surrounde d by circumste llar material that had been flash-ionized 



in the shock breakout (IFransson et al. l 1989). A few years after explosion, this 



material could be resolved as a triple-ring nebula consisting of an inner ring, 
located at a distance of 0.67 light-years and moving out with a velocity of 10 
km s^\ and two outer rings, located about three times further out and moving 
three times faster. The kinematics suggest a common origin in an ejection 
event some 20,000 years ago, the cause of which is still unclear. 

The shock, travelling at ~25,000 kms^^ reached the inner ring in 1995. 
A dramatic display of radiation, from X-rays to the radio regime, has since 
followed, caused both by the shock transmitted into the ring, and a reverse 
shock created by the impact. The emission has steadily increased over the 
years as the shock has penetrated into denser parts of the ring. In 2010, the 
X-ray curve flattened, p ossibly suggest ing that the densest parts of the ring 



have now been passed (Parket al.ll201lh . 



Recent calculations show that the X-ray emis sion from the ring started to 



influence the ejecta in 2001 (ILarsson et al.ll201 11) . Most of the X-rays hitting 



the ejecta are estimated to be absorbed in the hydrogen envelope, increasing 
the flux in hydrogen and helium lines. 

The ejecta itself, mostly confined to within 5000 kms^\ will start impact- 
ing the ring around 2015. A significant fraction of the 10^^ ergs of kinetic 
energy will then be converted to radiation (more than what has been emitted 
so far), however over a time-span of several decades. 
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3.6.4 Missing compact object 

A total of 24 neutrinos emitted in the collapse were detected by the 
Kamiokande II, 1MB and Baksan neutrino detectors, representing the only 
neutrinos detected from outside the solar system. The neutrino detection 
gave spectacular confirmation of the core-collapse model, with the number 
of neutrinos and their energy distribution being i n excellent ag reement 
with quantitative models for neutron star formation (IMcCrayl 120071) . That a 
neutron star was formed rather than a black hole is also suggested by the 
inferred explosion energy and the ^^Ni mass, which correspond to a rnass-cu t 
at ~ 1.6 Mq for any reasonable progenitor model ( Hashimoto et al. 19891) . 
well below the upper mass limit of 2 — 3 Mq for neutron stars. 

However, as of 201 1 no sign of the neutron star has been seen. It is possible 
that an optically thick dust cloud happens to obscure it. If the neutron star 
would be a pulsa r, i.e. magnetiz ed and spinning, it probably would have been 
detected by now dMcCraylboOTI) . 

Could a (non- spinnin g) neutron star in fluence the SN spectra? Neutron star 
cooling models by Tsuruta et all (l2002h show a total thermal luminosity of 
about 10 L0 for a decade-old neutron star. But the energy input by radioac- 
tivity is never below 500 Lq for the first 50 years, assuming a ^"^Ti mass of 
1 X 10^"* M0. Thus, the SN spectrum will not be noticeably influenced by 
the thermal emission of the neutron star. This is an important check for the 
validity of the '^^Ti-mass derived in paper II. 
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I now proceed to describe the main physical processes occurring in SNe, and 
how the state of the gas (temperature, ionization and excitation) can be cal- 
culated. The power source in the nebular phase is radioactivity, which heats, 
ionizes, and excites the gas, which then emits radiation. This radiation is partly 
reabsorbed, and partly escapes for us to observe. A schematic diagram show- 
ing the relevant processes is presented in Fig. 14.11 
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Figure 4.1: Schematic overview of the physical processes occurring in SNe. Radiative 
processes involving the internal UVOIR field are drawn as red, dashed arrows. 



4.1 Radioactive decay 

In the first step of a radioactive decay, a nuclear transmutation changes the 
number of protons in the nucleus. This occurs either by emission of an a- 
particle (a-decay), an electron plus an electron anti-neutrino (j8 "-decay), a 
positron plus an electron neutrino (j3+-decay), or by electron capture (EC), 
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where one of the orbiting electrons (usually from the K or L shell) is captured 
to convert one of the protons to a neutron, and the energy is usually emitted 
as an electron neutrino. The re-filling of the inner electron orbit occurs by 
emission of an X-ray or by the ejection of an electron (Auger electron). Apart 
from the decay processes mentioned above, some rare decays also involve the 
emission of a neutron or a proton. 

The new nucleus is usually formed in an excited configuration. It decays 
to the ground state, usually in several steps, by emission of y-rays (y-decay) 
and/or by internal conversions (IC), where the excitation energy is transferred 
to an orbiting electron, which is ejected. Following internal conversions, an 
X-ray or an Auger electron is also produced, just as after electron captures. 

A radioactive decay channel is thus associated with emission of 7-rays, X- 
rays, electrons, and positrons, all with individual energy distributions. Some 
channels may emit all of these particles, while others may emit only one or 
two types. Table |4~T] lists the fraction of the decay energy going to different 
particles for the most common radioactive isotopes in SNe. Note that ^^Co 
and ^^Sc are the decay products of ^^Ni and "^^Ti, respectively. One should 
observe that channels carrying only a small part of the total energy can still 
be important if the trapping of the decay product is more efficient than for the 
other ones. 



Isotope 


Gammas 


Positrons 


Electrons 


X-rays 




100% 





<0.1% 


<0.1% 


^^Co 


96.5% 


3.4% 


<0.1% 


<0.1% 


"Co 


86.4% 





11.3% 


1.8% 




92.7% 





6.7% 


0.49% 


^4Sc 


66.3% 


33.7% 


<0.1% 


<0.1% 



Table 4.1: Fraction of tiie decay energy (excluding the neutrinos) carried by various 
particles, for the most important radionuclides in SNe. The electron rest mass is not 
counted in the energy budg et, whereas t he posi t ron re s t mass is, being in cluded in the 
gamma column. Data from \Endl \jundA \l99^) . \Seiten7.ahk \201 7|) . 



4.2 Deposition of decay products 
4.2.1 Gamma-rays 

Gamma-rays see an opacity dominated by Compton scattering by free and 
bound electrons. Because the 7-ray energy (~ 1 MeV) is much higher than 
the binding energy of the electrons (< 1 keV), the 7-rays do not distinguish 
between bound and free electrons, and the deposition is therefore indepen- 
dent of the ionization/excitation state of the gas. The total cross section is 
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given by the Klein-Nishina formula, which gives the Thomson value {Ot = 
6.65 X lO^^^cm^) at low energies, but decreases monotonically with increas- 
ing energy. At 1 MeV the cross section is about 1/3 Gt, corresponding to an 
absorption coefficient 

Ky{l MeV) ~ ^Grm^'Y, = 0.066^ cm^g-i , (4.1) 

where Yg = Z/A is the number of electrons per nucleon. Most zones have 
Yg ~ 0.5, but the hydrogen zone has Y^ = 0.86, or lower if helium has been 
dredged up (Sect l2.3l ). 

For y-rays over 1.022 MeV, pair production can also occur in the vicinity 
of a nucleus. As the 7-rays are down-scattered below ~100 keV, photoelectric 
absorption adds to the opacity as well. 

Effective opacity 

Several authors (see below) have found that transport with a purely absorptive, 
effective (gray) opacity 

4^^-0.03^ cmV^ (4.2) 

gives an accurate approximation to the exact solution for the 7-ray transfer. 
The value is lower than the Klein-Nishina coefficient (Eq. 14.11 ) because this 
opacity represents pure absorption with no scatterings. The optimal value of 
Ky^^ has s ome dependen c y on t h e geome t ry of the pr oblem, but for typica l 



geometries IColgate et all (Il980l) . lAxelrodl (Il980l) . and IWooslev et a T (198ol) 



found values in the range k'/^ = 0.028 - 0.033fK /0.5). ISwartz et al.. (19951) 



made the important realization that the optimal value is time-dependent, find 

4" 



ing Ky-^^ = 0.033(ye/0.5) as long as multiple scatterings occur, falling to 



Ky^^ = 0.025 (7e/0.5) when zero or one scatterings occur. ISwartz et al.l(ll995h 



^7 

also provided a nice explanation why the effective opacity works so well 
for approximating the Compton scattering process; the interactions are such 
that forward scatterings remove little energy from the 7-ray, while strongly 
direction-changing ones remove a lot. Thus, the 7-ray travels along approxi- 
mately straight paths with little energy loss until a strong interaction occurs, 
at which point it loses most of its energy. This situation is obviously well 
approximated by a purely absorptive transfer. 

In this thesis, gray transport with Eq. l4.2l is used. One should be aware that 
this treatment is only accurate to 10 — 20%. However, as we in general do not 
know the detailed distribution of the ejecta zones within the SN, performing 
exact Compton calculations is not very meaningful. 

A homogeneous sphere of mass M and expansion velocity V reaches optical 
depth unity to the gamma-rays at a time 

C = 460 days ( ^'^ ( , ^ ' . (4.3) 
^ ^ \lOMcJ V3000km/sy ^ 
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Since the expansion velocity V can be estimated from the line widths, the 
epoch when the bolometric light-curve starts decreasing faster than g^'/'^seco 
can therefore give us an estimate of the ejecta mass M. Fig. 14.21 shows the 
bolometric light-curve of SN 1987A, normalized to the decay by 0.07 of 
^^Co. We see that the trapping reaches 1 — e^' = 0.63 at about 530 days after 
explosion. From Eq. 14.31 with W = 3000 km s^\ the ejecta mass can then be 
estimated as 13 M©, in good agreement with estimates from more detailed 
modeling. The assumption of constant density, as well as which V to use, are 
the limitations of the method, but Eq. I4.3l can still be used to tell us the ejecta 
mass to within a factor two or so. 
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Figure 4.2: The bolometric ligh t-curve of SN 1987 A norm alized to the powering 
by 0.07 Mq of ^''Co. Data from ISuntzeff & Bouche3 ( Il99d) . Time is in days since 
explosion. 



The 7-ray deposition as function of radius in a uniform sphere is plotted in 
Fig. l4.3l for a few optical depths. The difference between the outer edge and the 
center never exceeds a factor of two. From the perspective of y-ray deposition, 
it is therefore reasonable to use single-component zones in the core where the 
^^Ni is distributed. Such models are used throughout this thesis. 



4.2.2 Electrons and positrons 

The electrons and positrons emitted in the radioactive decays lose their energy 
by ionizing and exciting the gas. Since their energy is much higher than the 
excitation and ionization potentials, one may use the Bethe approximation for 
both ionization and excitation cross sections (lBethelll930h . The cross section 
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Figure 4.3: Gamma-ray deposition (arbi trary units) as function of r adius in a uniform 
sphere for different optical depths. From lKozma & Franssonl(ll992[ KF92 hereafter). 



for ion /, transition j, is 



Ocoujj iE) = -^ {Aij lnE + Bij) E I , 



(4.4) 



where E is the energy of the colliding electron, ao is the Bohr radius, and A/y 
and Bij are constants, depending only on level energies and A-values. For 
lower energies, experimental data for the collision cross sections are needed. 
KF92 provides a summary of the atomic data we use. A,y is, for bound-bound 
transitions, proportional to the oscillator strength fij. Allowed transitions 
(/;;■ ~ 1) will therefore in general be more important than forbidden ones 
{fij <C 1). Most elements have a continuum oscillator strength that is similar 
to the sum of excitation oscillator strengths, so about as much energy goes 
into ionizing as exciting the gas by these primary particles. However, the 
total ionization and excitation rates are influenced mainly by the secondary 
electrons ejected in each ionization, described more in Sect. 14.31 

Since 1 MeV ~ 10^ times the ionization/excitation potentials (~ 10 eV), 
some 10^ scatterings are involved before the kinetic energy has been converted 
to ionization and excitation energy. The time-scale for the down- scattering 
process can be shown to be, for a pure hydrogen nebula, (■Axelrod.1980,) 



t slowdown ~ 10 days 



M 



10 Mr. 



V 



3000 km s 



-1 



100 days 



(4.5) 



which stays short for many decades. 

An effective, gray opacity c a n be u sed f or the lep t on tr ansport just as for the 
/-ray transport. IColgate et al.l (119801) and lAxelrodl (119801) found good values 
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for the absorption coefficient to be 

KfJ = {l-lO)^ cmV\ (4.6) 



from which it is clear that the electrons and positrons can travel only ~ 1 /300 
the distance of the y-rays. A SN (without a magnetic field) becomes optically 
thin to the leptons at a time 

which shows that they remain fully trapped for several decades, and further- 
more we can assume on-the-spot absorption for the first year or so. A disor- 
dered magnetic field will further increase the trapping by locking the particles 
in La rmor orbits. Th e cyclotron radius Rc, relative to the radius of the nebula 
R, is (lAxelrodlll980l) 



R \10-('GJ V3000kms-V VlOOd/ ^ ^ 

A magnetic field even at the interstellar level (B ~ 10^^ G) will therefore 
produce local trapping of the positrons, and the trapping will increase with 
time. In paper II, we find that such a field appears to be present in the ejecta 
of SN 1987A. 



4.2.3 X-rays 

The X-ray opacity is dominated by photoelectric absorption on inner-shell 
electrons of metals such as C, O, Si, and Fe. The total photoionization 
cross section at 1 keV , for s olar abundances, is 3 x 10^^^ cm^ (e.g. 

19831) . The trapping time-scale can then be 



Morrison & McCammon 
computed as 



trap 
'■X-rays 



540 years 



M 



10 



1/2 



V 

3000 km s"' 



-1 



Ex 
1 keV 



-3/2 



(4.9) 



Apart from radioactivity, the ejecta is exposed to (at least) two other sources 
of X-rays; the new-bom neutron star and the shock-heated circumstellar 
gas. The neutron s tar will have a temperature of a few million degrees (e.g. 
Tsuruta et al ] l2002h . emitting photons with an average energy of 3.82 kgT ~1 
keV. However, as we saw in Sect. 13.6.41 this flux is too weak to influence the 
SN spectr um. X-ray input from circumstellar gas may be more important. 
Recently, iLarsson et al.l (1201 ih showed that the X-rays from the forward 



shock in SN 1987A are influencing the ejecta conditions at late times. 
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4.2.4 Deposition time-evolution 

Fig. 14.41 shows the fraction of the total energy deposited in the H, He, Fe/He, 
and 0/Ne/Mg zones as function of time for the 12 M0 model in paper IV. 
The optical depth to the 7-rays is initially quite high in the Fe/He-clumps, but 
as the trapping in these clumps decreases, so does the fraction deposited here. 
After 300-400 days, the deposition in the Fe/He clumps becomes dominated 
by positrons. Since the 7-ray trapping in the other core zones decreases as 
t~^, their curves continue decreasing while the Fe/He curve turns up. The 
deposition in the H-zone is a monotonically increasing function with time as 
long as the total 7-ray deposition exceeds the total positron deposition. This is 
because the deposition in the H-envelope falls off more slowly than f^^, as a 
larger and larger fraction of the 7-rays escape from the core. Clearly, the best 
time to study the Fe/He (and Si/S) zone is for times later than ~ 1000 days, 
something exploited in Papers I and II. The best time to study the oxygen 
zones is instead as early as possible in the nebular phase, the basis of papers 
III and IV. 




100 200 300 400 500 600 700 

Time [days] 

Figure 4.4: Fraction of the total deposition going to the H, He, Fe/He, and O/Ne/Mg 
zones as function of time, for the 12 M© model in paper IV. 
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4.3 Ionizations, excitations, and heating produced by 
radioactivity 

Each ionization by a y-ray or high-energy lepton leads to the creation of a 
new high-energy electron. These secondary electrons will excite and ionize 
the gas just as the primary ones. Furthermore, once the electrons are slowed 
down to keV energies, interaction with the thermal electron pool becomes 
important, leading to heating. In this way, the radioactive decays lead to ion- 
izations, excitations, and heating of the gas. At each point in time and space 
we will have a distribution describing the number of non-thermal electrons per 
energy interval. If this distribution can be determined, the rates for all these 
processes can be calculated by integration over the relevant cross sections. Of 
course, these rates are precisely what determines the distribution, so the sys- 
ter n has to b e solv ed self-consistently. It was first solved in the SN context 
by lAxelrodl (Il980h . using the Continuous Slowing Down approximation (in 
which the electrons are assurt ied to lose o n ly a small part of the ir energy in 
each interaction), and later bv lLucyl(ll99ll) . IXu & McCrayl (Il99lh and KF92, 
using exact solutions. These authors found that the solution is insensitive to 
the exact properties (type and energy) of the original input particles, and to the 
density. The only major determinants are the chemical composition (the ion 
abundances) and the electron fraction. Fig. 14.51 shows the solutions for pure 
oxygen and iron plasmas. 
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Figure 4.5: The fraction of radioactive powering going into heating, ionization, and 
excitation in a pure oxygen plasma (left) and a pure iron plasma (right), as function of 
the electron fraction Xg. From KF92. 



Throughout this thesis, the KF92 code (with some updated atomic data) has 
been used for computing the non-thermal ionization, excitation, and heating 
rates. To save computation time, only ions with number fractions over 10^^ 
were included in the solutions. For the others, the average ionization rate for 
the corresponding ionization degree was used, whereas the non-thermal exci- 
tation rates were put to zero. 



4.4 Temperature 
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4.4 Temperature 

Many of the processes important for the thermal structure of S Ne are the same 



as in H I regions, discussed by iDalgamo & McCrayi (1 19721) . Both environ 



ments are low-ionization nebulae heated mainly by non-thermal processes. 

The first law of thermodynamics states that the energy content ?7 of a parcel 
of gas changes as 

^ = G(0-MO^, (4.10) 

where Q{t) is the net heating rate of the parcel, Q{t) = H{t) — C{t), where 
H{t) is the heating rate and C{t) is the cooling rate, p{t) is the pressure, and 
V{t) is the volume. Cooling may, in general, occur by conduction, convec- 
tion, or radiation, but in SNe only radiative cooling is important. For an ideal 
mono-atomic gas we have U (t) = 3/2N{t)kBT{t), v/hereN{t) is the total num- 
ber of particles (atoms and electrons) at time t. If Na is the total number of 
atoms, then N{t) =Na{l +Xe(t)), where Xe{t) is the electron fraction. Further, 
p{t) = n{t)kBT{t), where n{t) =N{t) /'V{t). We also have, for a homologously 
expanding gas, V{t) = Vo(?Ao)"^ so dV{t)/dt = 3V{t)/t. Eq. l4TT0l can then 
be written as 

dT{t) h{t)-c{t) 2T{t) T{t) dx,{t) 



dt Usnit) t l+Xe{t) dt 



(4.11) 



where h{t) and c{t) are now the heating and cooling rates per unit volume. The 
last term is usually much smaller than the other ones and may be neglected. 

The cooling c{t) is usually dominated by coUisional excitations giving 
atomic line radiation, just as in H I and H II regions. This line cooling 
increases with temperature, since a higher temperature means that a larger 
fraction of the electrons have high enough energy to excite any given 
transition. If coUisional deexcitations are unimportant, c{t) also increases 
with density as nc{t)ne{t), where nc{t) is the number density of cooling atoms 
and is the electron number density. 

Since also the adiabatic cooling term, —2T{t)/t, increases with tempera- 
ture, a cold gas has weak cooling compared to the (presumed temperature- 
independent) heating. If heating of this gas is initiated, only the h{t) term 
in Eg. 14. Ill will be important and the temperature will increase. For each in- 
crease, c{t) will increase, and so h{t)—c{t) will become smaller. The adiabatic 
cooling increases as well if T grows more quickly than oc t. Eventually, heating 
and cooling will match, making dT /dt small, and the temperature stabilizes. 

The thermal equilibrium approximation corresponds to setting 

h{t)=c{t), (4.12) 

turning the differential equation 14.111 into an algebraic equation. Since both 
heating and cooling are radiative (if we count radioactive decay particles as ra- 
diation), this condition corresponds to radiative equilibrium (in the co-moving 



36 



Physical processes 



frame). From Eq. 14.111 the approximation is good if (neglecting the last term) 

Ht) c(t) 2T(t) 
^ ^ ^ ^ > T{t) + — ^, (4.13) 



jkBn{t) jkBn{t) 
which can be written as 

fheat(0,fcool(0 < fchange(0 , (4.14) 

where fchange(0 is the time-scale for change in physical conditions (density, 
temperature) 

frit') 2\^^ 

?change(0=(^^ + 7j , (4-15) 

and fheat(0 and fcooi(0 are the time-scales for adding and removing the thermal 
energy content, respectively, at a given time: 

lksn{t)T{t) \ksT{t) 
\kBn{t)T{t) \kBT{t) 

'-'^'^ = c{t) = -cJtT ' ^^-''^ 

where Hpp{t) and Cpp(f) are the heating and cooling rates per particle. By 
inspection of Eg. 14. Ill it is is clear that either both conditions in Eq. 14. 141 will 
be fulfilled, or none of them will. It is therefore enough to compare one of the 
time- scales, usually chosen as fcooi, to ^change- 

The temperature will be able to significantly change no faster than on the 
shorter of the radioactive decay time- scale and the expansion time- scale. If 
the latter is the shorter one, we can approximate j/r ~ T{t)/t, which gives 
tchange{t) = t /3. If the radioactivc time-scale Xdecay is shorter, we can write 
T{t)/T{t) ~ l/tdecay, and then tchange ~ Trfecay- Thermal equiUbrium corre- 
sponds to to having the cooling time-scale (or equivalently, the heating time- 
scale) shorter than the smallest value of f/3 and Tdecay 

The code uses the thermal equilibrium approximation to compute the tem- 
perature at any given point and time. As Fig. 14.61 illustrates with an example 
calculation, the steady-state approximation can cause both under and overes- 
timates of the actual temperature. 



4.4.1 Heating rates 

Heating may in general occur by radioactivity (non-thermal electron colli- 
sions), photoionizations, line absorptions followed by coUisional deexcita- 
tions, free-free absorptions, exothermic charge transfer reactions, and Penning 
ionizations: 

h = hnt + hpi + hcd+hff + hct + hpe . (4.18) 
These calculation of these are described in paper II (Sect. 2.3.3). 
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Figure 4.6: Example calculation of the temperature evolution T{t) for an exact so- 
lution (solid line), and a thermal equilibrium approximation (dashed line). The ap- 
proximate solution starts breaking down at ^ 200 days, giving an overestimate of the 
temperature (and thereby the thermal emission) from ^ 200-500 days. The loss of 
cooling efficiency (in the exact solution) means that deposited energy is stored up in 
the gas. This causes the approximate solution to underestimate the temperature (and 
thereby the thermal emission) at late times. 

4.4.2 Cooling rates 

Cooling may in general occur by atomic line cooling, recombinations, free- 
free emission, endothermic charge transfer reactions, molecular line cooling, 
and dust cooling 

C = Cline + Crec + Cff + C^t + C,nol + Cdust ■ (4- 19) 

These are also described in paper II, but I make some additional comments on 
them here. 

Line cooling cum 

The cooling in a particular line is in general not equivalent to the emission in 
the line, as population of the upper level may occur also by other processes 
than by energy transfer from the thermal electron pool. Instead we have to 
consider how the thermal electrons lose and gain energy by exciting and de- 
exciting the transition. Letting / and u denote the lower and upper level: 

c|L = A£/„«. [mqiu (T) - n^qui {T)] , (4.20) 

where AE/^ is the transition energy, and qiu{T) and qui{T) are the rates for 
upward and downward collisions, related by 

qiu{T) = ^-^qu,{T)e-'^'""^^\ (4.21) 
gl 
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where gu and gi are the statistical weights. We can therefore write 



Ju 
''line 



^^lunequlij) 



Si 



(4.22) 



If the upper level population ?i„ exceeds ^me^^'"!^'^'^ (its LTE value «* 
in the two-level approximation), the transition therefore induces net heating 
instead of cooling. This is often the case for low-lying H I and He I transitions, 
where non-thermal excitat ions and recombinations over-populate the excited 
states dde Kool et al.lll998l) . 

The coUisional cro ss sections typically have a 1/V/ dependence 

( Osterbrock & Ferland 2006h . where Ve is the velocity of the thermal 
electrons. The effect is that, even though more encounters occur per second in 
a hotter gas (in proportion to Vg), there are actually fewer collisions occurring 
(for a given electron energy) by a factor Vg °<- T^l"^. The qiu coefficients 
therefore have a temperature-dependence 



qiu{T)=T-'I^Ciu{T) 



(4.23) 



where C/„(r) are transitions-specific functions, weakly dependent on T . 

Now, consider a two-level atom with only collisions and radiative decays 
occurring. Clearly, the cooling in this particular situation is equal to the emis- 
sivity in the line (since no other process is allowed to populate the upper level): 



Ju 
" line 



riuAEiuAuip, 



(4.24) 



where j8„/ is the escape probability. The upper level population has its max- 
imum value in LTE: 

(4.25) 



n„* 



gl 



But if we insert this into Eq. 14.221 we get c|.^^ = 0! The resolution to this 
paradox lies in the realization that LTE and any net radiative losses are not 
compatible; strict LTE requires that all deexcitations are coUisional. The ra- 
diative cooling corresponds, equation-wise, to the potentially small deviation 
of Hu from its LTE value («„*) 



1 



AulPul 



Ju 
" line 



AEiuUequiit) nl 



'ul 



'qul{T)n, 



quine 
: nlAEi^AuiPui . 



(4.26) 
(4.27) 



The ratio Aui^ui / quing may become very small, and the actual computation 
of the cooling rate then constitutes taking the difference between two much 
large r numbers. This may potentially cause problems by numerical cancella- 
tion (lHeathlll997l) . The machine precision at double precision is ~ 10^^. The 
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minimum value of A^iPui /luine is of order 1 /ne, so for > 10^^ cm'', numeri- 
cal cooling calculations by this method break down, and one instead has to use 
the emissivity combined with some method to estimate the non-thermal con- 
tribution. This situation, fortunately, occurs well above the electron number 
densities in nebular phase SN environments. Potential problems may instead 
arise by insufficient accuracy in the NLTE levels computations. If we com- 
pute these to an accuracy of 1%, but the cooling corresponds to a difference of 
0.1% between upward and downward collision rates, our cooling rates become 
wrong! This situation is usually salvaged by the fact that the low-lying levels 
that are important for the cooling converge more rapidly than the high-lying 
levels. When global convergence is reached at the 1% level, the accuracy in 
the important levels is therefore much higher. 

An important point is that line cooling cannot be quenched by increasing 
the density so that the coUisional deexcitations increase; this effect is always 
superseded by an increase in the number of upward collisions. The maximum 
cooling rates are achieved in LTE. 

Recombination cooling Crec 

The recombination cooling rate per unit volume equals the recombination rate 
times the average energy of the recombining electrons. This value is ~ Q.^ksT 



for hydrogen at 5000 K (lOsterbrock & Ferlandll2006h . which we use for all 
atoms. Thus, 

Crec = neY,OCkiT)nkOMBT , (4.28) 

k 

where ak{T) is the total recombination rate for ion k. 

Free-free cooling Cff 

The free-free cooling is given by (e.g. lOsterbrock & Ferlandll2006r) 



Cff = 1.42 X 10-27„^ri/2£z2g|«k , (4.29) 



where summation occurs over all ions k. We follow Ide Kool et al.l (Il998h and 
use a constant Gaunt factor g'^ = 1-3. 

Charge transfer cooling Cct 

Charge transfer cooling may occur by endothermic reactions, which are in- 
cluded as negative terms in the net charge transfer heating rates (Sec. l4.4.TI ). 

Molecular cooling Cmoi 

Molecules are efficient coolers due to their rich level structure. The models 
in this thesis do not include calculation of the molecular abundances, and so 
we cannot compute the molecular cooling rates explicitly. In papers I and II, 
where the spectra anyway are rather insensitive to the temperature, we set 
Cmoi = 0. In papers III and IV, we attempt a more realistic treatment be setting 
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Cmoi = ^ in the O/Si/S and O/C zones, where SiO and CO are expected to 
form. 

Significant work has already been done for the inclusion of molecules in 
the code, and will hopefully be achieved in the near future. 

Dust cooling c^ust 

Dust can be a coolant by inelastic colhsions of the electrons with the grain 
surfaces. We do not include this process in the code, so Cdust = 0. 



4.5 Ionization 

In nebular phase SNe, the gas is far from LTE, so we have to compute the 
ionization balance from the rate equations. Letting xjj denote the fraction of 
element j in ionization state /, the time-derivative of xj^i is 

= r^-,,-i^7,i-i - {¥j,i + + x(fj,i+iXj,i+i , (4.30) 

where Fjj is the ionization rate per particle and is the total recombina- 
tion rate per particle. The ionization equilibrium approximation corresponds 
to setting the time-derivative to zero, which gives equations of the form 

Tj^i-iXj^i-i = Yj,iXj,i . (4.31) 
Together with the number conservation equations (one for each element) 

£x,y = l, (4.32) 

(=0 

and the charge conservation equation 

z(i) 

"e = II"yH'^7V' (4-33) 

(=0 

the ionization balance, i.e. the solutions for all xjj, can be computed. The sys- 
tem of equations is non-linear because the ionization and recombination rates 
generally depend on the ion abundances themselves (through dependence on 
ion densities for charge transfer reactions and on the electron number density 
for recombination rates), and is solved by the Newton-Raphson method. 
The ionization equilibrium approximation is valid for 

<■ rjj_ixjj_i + xifj^i+ixj,i+i 

(4.34) 
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with the same reasoning as in Sect. 14.41 either both conditions will be fulfilled, 
or none of them will. We can therefore choose the one easier to test. This is, 
usually, the second one. We can then write 

fo'ut(0 <■ (^^^^) ~ min[f, Tdecay] , (4.35) 



where 



Since t^L < f/ec, where 



t) = (4.36) 



tgit) = (4.37) 
rj,i 

(4.38) 

ionization equilibrium is fulfilled if tree ^ min[f, Xdecay]- Assuming the recom- 
bination rate to be dominated by electron recombinations, we have 

fj' = — i — (4 39) 

The recombination rate is of order ajj ~ 10^^^ cm^s^' at a few thousand 
degrees, so the recombination time-scale is, for Xe = 10^^ 

M V V / Xe 



The core region of SNe will therefore have a recombination time scale shorter 
than t and idecay for a decade or so. The density gradient in the outer envelope 
is so steep, however, tha t the recombination time- scale there becomes long 
already after 1 — 2 years dpransson & Kozma 1993 ). 

The ionization rates will decline exponentially with the radioactivity, 
whereas the recombination rates will evolve roughly with density as For 
small t, the t^^ decline is more rapid, but at later times the exponential decay 
is the more rapid. The result is that the ionization degree reaches a maximum 
after some time, and then declines. If we approximate TjjXjfi = Ce^'^''^''"-'"' 
and ajj = constant, the solution for the ionization fraction in a single-species 
gas is 

Xe{t) = Cf'-f/2^''--Vf3/2 ^ (441) 

which reaches its maximum at t = 'iXdecay, so we can expect the ionization 
degree to peak around ~ 1 year (3-111 days), independent of the ^^Ni mass 
as well as the ejecta mass and velocity. 

Fig. 14.71 shows an example calculation using the exact equation (Eq. 14.301) 
and the ionization equilibrium approximation (Eq. 14.311 ). Just as for the tem- 
perature, the ionization fraction is initially overestimated, but then underesti- 
mated. 
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Figure 4.7: Example calculation of the ionization fraction of a two-level system using 
the exact solution (solid blue), and the ionization equilibrium approximation (dashed 
red). 



4.5.1 Ionization rates 

Ionizations can occur by radioactivity (non-thermal electron collisions), pho- 
toionizations, or by charge transfer: 

Fjj = rf,+rf.+rfj. (4.42) 

In SN nebular phase environments, the temperature and density are too low 
for coUisional ionizations by thermal electrons to be important. 

Non-thermal ionization rates (T^^i) 

These are obtained from the solution to the degradation of radioactive prod- 
ucts described in Sect. 14.31 In this model, ionizations may occur for any of the 
outer shell electrons, as well as for electrons in the shell inside that. The total 
ionization rate should be the sum over the individual contribution from these. 
Letting Ejj denote the energy per second and unit volume going into ioniza- 
tion of ion j,i, and letting nj" denote the ground-multiplet number density, 
and Xji the ionization potential of bound electron k, the ionization rate for the 
ground state as a whole is 

pnt,exact _ j,i j,i _ .r,. 



where 



We currently approximate the averaged ionization potential Xjj with the out- 
ermost electron ground state potential, Xjj = and leave the ionized 
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atoms in their ground states. The assumption of only ground multiplet ion- 
izations is a good approximation in the nebular phase when few atoms are in 
excited states. We thus put 

F ■ 

r"t. = II (4 45) 

7,' gm ^valence ' V^-^-^y 



An exact treatment would involve using the exact expression ("Eg. 14.441 ). but 
then also creating some ions with electron gaps, and following the subsequent 
fluorescence and Auger events. 

Photoionization rates fr^',) 

The photoionization rates are, in general, related to the mean intensity of the 
radiation field /y (see Chapter 5). If a fraction yi of a particular ion are in 
excitation state /, the contribution by level / to the photoionization rate of the 
ion is 

rf^, =y,r o,y^dv, (4.46) 
■'" JvO{l) hv 

where a/ y is the photoionization cross section, and Vo(/) is the threshold. The 
total ionization rate is thus 

^i = Ly'r ^^^'^^^^ (4.47) 

One major advancement in this thesis is the development of a radiative transfer 
calculation to be able to compute these rates, which is done by direct counting 
of the photoionization events in the Monte Carlo simulations. Details of these 
computations can be found in paper II (Sect. 2.3.1) and paper III (Appendix 
A). The photoionization cross sections we use are the ones for ionization to 
the ground state in the ionized atom. One exception is for O I ionizations, 
where ionizations to the first few excited states are included. 

Charge transfer rates (^\) 

Charge transfer (CT) refers to reactions where an electron is transferred from 
one atom/ion to another. Letting superscripts denote the level of ionization, 
ionization of X' by F' can be written 

X' + F'" X'+i + F''"i (4.48) 

The total ionization rate is the sum of the rates with each individual 
atom/ion. 

f,i' 

The reaction rates can be large, ^f-'T-, ~ 10^^ cm-'s^^ especially if, in addi- 
tion to fulfillment of selection rules, the energy defect of the reaction is small. 
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< 0.6 eV (lMeliuslll974l) . Often, this is fulfilled by reactions occurring to some 



excited state in one of the reaction products. 

Both the non-thermal and photoelectric ionization rates can be expected to 
decline exponentially with the radioactive decay. Charge transfer ionization 
rates, however, are proportional to the ion density, which Eg. 14.4 II combined 
with n oc t^^ shows to decline as ~ g^'/^^f^^/^. Thus, we have roughly 

. oc e'/^'r'^^ , (4.50) 

which has a minimum att = 3t. The implication is that charge transfer has its 
smallest impact after about one year, but plays an increasingly important role 
for the ionization balance at late times. 

In highly ionized plasmas, charge transfer is less important because of the 
Coulomb repulsion between ions. But in the low-ionization environment of 
nebular phase SNe, they are critical, and lack of atomic data is a major obsta- 
cle. 



4.5.2 Recombination rates 

Recombinations can occur by radiative recombinations, dielectronic recombi- 
nations, or charge transfer. 

VA,-,,- = vaJ:,- + V^HV^} (4.51) 

Collisional recombinations are unimportant. Stimulated recombinations can 
be treated as negative photoionizations, but have not been considered in this 
thesis. 

Radiative recombinations Yji 
The radiative recombination rate is 

xifff'^' = ajj{T)n, , (4.52) 

where CCjj{T) is the radiative recombination coefficient for the relevant ion, 
averaged over the internal states of the ion: 

^j,iiT)='£yiaj.iAT) (4.53) 

Since most atoms will be in the ground multiplet, we approximate cCjj{T) = 

\irli = aj,i,g,n{T)ne (4.54) 
The temperature dependence for the recombination coefficient is typically 
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Dielectronic recombinations \f/j] 

Dielectronic recombinations involve energy transfer to one of the bound elec- 
trons rather than emission of a free-bound photon; it is thus radiation-less. 
The atom is left in a state with two excited electrons. One of them may be 
ejected again in a so called autoionization, or they may both cascade to the 
ground state. Dielectronic recombination rates are usually significant only at 
very high temperatures (~ 10^ K), but can be compar able to the radiative re 



combination rates for some ions at low temperature (iNussbaumer & Storey 



1983h . They are currently not included in the code: 

Wfj = 0. (4.55) 

Charge transfer recombinations 

The charge transfer recombination rate is 

As mentioned above, the reaction rates may be large, ^'^^(f^ ~ 10^^ cm^s^\ 
compared to ay , ~ 10^ '•^ cm^ s^' for recombinations with free electrons. In 
a neutral gas, where the number density of neutral species is much larger than 
it is clear that v^' may be orders of magnitudes larger than . Charge 
transfer recombination therefore dominates the recombination rate for many 
ions, especially the ones with high ionization potential that can ionize many 
other elements. One example of this is oxygen ions in SN cores; large abun- 
dances of C I, Na I, Mg I, Si I, and others lead to charge transfer recom- 
binations dominating radiative ones, quenching oxygen recombination lines 
(paper 11). 



4.6 Excitation 

Just as for the ions, we can set up equations for the population of excited states 
within each ion. Letting denote the fraction of some element in ionization 
state / and excitation state /, we have (if someone knows how to write fun 
equations, let me know!) 

V=jLl I' I' / 

where the terms represent ionizations into the level (r,_i // , /), recombinations 
into the level (V^z+i,/',,,/), radiative and coUisional transitions from all other 
(internal) levels into the level (P/' /), radiative and coUisional rates out of the 
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level to all other (internal) levels (P/ ;')> ionizations from the level {Pij,i+ij'), 
and recombinations from the level (Yi,i,i-ij')- 

The excitation equilibrium approximation corresponds to setting the time- 
derivative to zero. Analogously with the thermal and ionization equilibrium, 
the approximation is valid if the deexcitation time-scale is short compared to 
t and Tdecay Sincc deexcitation rates are never smaller than ~ 10^^ s^' for 
excited states, this condition is always fulfilled for them. Thus, for given ion 
abundances, we can always compute the excitation structure in steady-state. 

The resulting system of equations would be linear (for a given n^) if all 
lines would be optically thin. In general, the system is non-linear because 
the deexcitation rates P/ // depend on the lower level populations for optically 
thick lines. We solve the system of equations by straight iteration; guess the 
solutions /, compute the deexcitation rates Piji, make a new solution for /, 
and iterative until convergence. 

To enhance the convergence stability, the model atoms should be kept as 
small as possible. Also reasonable execution times is a strong reason to keep 
the model atoms limited in size. All in aU, the construction of the model atoms 
is a complex optimisation problem, where a compromise between complete- 
ness, accuracy, stability, and efficiency has to be found. 



4.6.1 Transitions rates 

Photoionization and recombination rates (F,!//^) have been described earlier 
in the chapter. Here, I describe the internal transition rates (P; /')- Letting 
/ and u denote the lower and upper levels, the internal transition processes 
are spontaneous emission (Rf""'), stimulated emission (R^Jf), photoabsorp- 
tion (Rfj), thermal collisions (c'^ymai ^qhermai^^ non-thermal collisions 

//^non— thermal /-^non—thermal\. 



D Tispom , Tistim I /"'thermal , /-^mm— thermal /a cq\ 

Pud = Ru,l +^u,l +<-u,l +^u,l <.4.5«) 

D nabs I /-^thermal , r^mm— thermal /a cn\ 

Plu = Ri,u +^i,u +(-i,u ■ <^4.5y) 



Spontaneous emission rate R^^ 

The spontaneous emission rate is given by the Einstein A-coefficient. How- 
ever, in the Sobolev approximation, one cancels the contribution by this spon- 
taneous emission to the R"''^ and rates (see Chapter 5, Sect. 15.2.21 ). and 
instead uses 

K!:r=Auil5!n (4.60) 

where jS^^^ is the Sobolev escape probability (see Chapter 5). For some lines, 
P^i is replaced by an effective escape probability jS,^/^ , which includes absorp- 
tion in the continuum and/or in other lines. See Sect. l5.2.2l for more details. 
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Radiative absorption rates R"'^J 

The general expression for the radiative absorption rate is 

poo 

Rf^'''"" =B,,J^ Jy^ydV=B,J,, (4.61) 

As mentioned above, contributions by spontaneous emission in the line itself 
should be excluded in //„ in the Sobolev approximation. One may then show 
(see Sect. 15.2.21 ) that the expression to use is 

R^^=BiM (4.62) 

where Jy represent the radiation field in the blue wing of the line (where no 
Une photons exist). Our method for computing is described in paper III 
(Appendix A). 



Stimulated emission rate 



We obtain the stimulated rates from R'}"^ according to 

= ^Rl^^ (4.63) 
8u 

Thermal collision rates C[^^"""^ 

No selection rules apply for coUisional transitions, so they may be important 
for all lines. Typically only electron collisions are considered, as the thermal 
velocities, and hence the collision frequencies, are much lower for atoms and 
ions 

Vafnms ( nty \ 1 , . ^ a\ 

(4.64) 



Velecti-ons \Amp J 43 A^/^ 

so the colUsion rates scale as 

Cf™' (atoms) 



CJ^J^™"' (electrons) 43Ai/2 



(4.65) 



Collisions with atoms can thus be ignored in astrophysical plasmas of high 
enough ionization degree {Xe > 1/(43 A^/^)). In addition, collisions with 
atoms are less efficient than collisions with electrons and ions due to the 
long-range electric force (and ions are obviously no more common than 
electrons). Still, when falls below ~1%, hydrogen atom collisions (A = 1) 
may become important. 

The thermal collision rates are 

(^thermal ^ 8.629 X IQ-^n^r' '/^ , (4.66) 

§u 

^thermal ^thermal Si ^—MLui /hgT (4 67) 
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where T/ „ are the velocity-averaged collision strengths, which are slowly 
varying functions of T. When they are unknown, approximations have to be 
used. For allowed transitions, van Regemorter's formula is co mmonly used . 



For forbidden transitions, rates are often computed according to lAIlenI (119731) . 
where an effective collision strength of 1 is assumed. 

Non-thermal collision rates cfjf 

For the computation of the non-thermal excitation rates in KF92, the indi- 
vidual levels in the ground multiplets are generally assumed to be populated 
according to their statistical weights, which is true if IcbT » AE, where AE is 
the internal energy scale of the multiplet. For typical AE of ~0. 1 eV, this cor- 
responds to r » 400 K. For the very late-time models considered in papers I 
and II, this condition is not fulfilled. We therefore modified all cross sections 
by factors e^^*/*^^^ in the Spencer-Fano routine, where Ei is the excitation en- 
ergy of level /, which approximates the deviations of the level populations 
from the statistical weights limit. The result is more energy going into chan- 
nels from the low-lying states in the ground multiplet. 



4.7 Global convergence 



AxelrodI (119801) showed that the equation sets 14.12114371 and 14.571 (with the 
time-derivative set to zero) may be solved for global convergence by itera- 
tion. All systems are non-linear. The most commonly used solutio n technique 



for a set of non-linear equations is the Newton-Raphson method (iPress et al. 



1992D . Here we make a starting guess for the solution, linearize all equations 
at that point, and solve the resulting linear system. In other words, we replace 
fiix,y,..) = with fl'"{x,y, ..) = 0, where 

fl'^ix,y, ..) = Mxo,yo, ■■) + ^{x-xo) + ^iy-yo) + - (4-68) 

We then repeat the procedure for the new solution point. For this scheme to 
work, the linearized equations must have a solution at least in the same direc- 
tion at the real one. This is true in the vicinity of the true solution, but may not 
be true far away. The success of the Newton-Raphson method therefore relies 
on making good enough starting guesses. We use LTE solutions for ionization 
and excitation as starting guesses in the NLTE computations. 



4.7.1 Solving a set of linear equations 

The computation of each Newton-Raphson step requires the solution of a 
linear system of equations. The techniques for this are well known, and 
numerical packages are readily available for the task. But let us have a look 
at the estimated requirements of disk space and computation times for the 
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problem. 

A linear problem (Ax = b) with 10^ variables has a matrix A with 
10^ entries. Even storing that many numbers takes ~ 8 • 10^ ~ 1 GigaBytes. 
The number of floating point operations needed to solve a syste m of size 
using the standard method of LU decompositiorQ is ( Heath 1997 ) 



Nfp ~ ^A^^ , (4.69) 

which tells us that it is very expensive to include more atoms in the ionization 
balance and more levels in the excitation structure. It also tells us that, in 
general, it is a good idea to split the problem into chunks which are solved 
iteratively, rather than attempting to solve everything in one big calculation. 
There is however, a problem in that we intuitively can expect convergence 
difficulties to increase with increasing fragmentation of the equation set. We 
therefore need to strike a balance in the amount of splitting that is done. 

For an atom size of N = 10^ we have Nfp = 3- 10^. On a 1 GHz CPU, the 
computation time is then ~ 0.03 second^. Since we need to solve for ~ 10^ 
atoms in ~ 10 zones over perhaps ~ 10^ iterations, the total CPU requirements 
for the statistical equilibrium is then ~ 1 hour. We parallelize the computations 
over the zones, reducing the run times to ~ 0.1 hours. 



4.7.2 Convergence 

The stability of a Newton-Raphson algorithm often increases if we use a 
damped step-size rather than the full one. This means that the step imple- 
mented is some fraction h of the one computed from the solution to f!'" = 0. 
We then reduce the risk of over- sho oting t he solu tion to a non-convergent 
regime. In t he hterature, w e find that ILucyi (12003 ah uses a damped step size 
h = 0.8, and lLucvl (Il999bh uses h = 0.5 — 0.8, furthermore with a temperature 
step size limit of AT = 200 K. The convergence can be further enhanced by 
making h and AT adjustable to smaller values if the step sizes are not decreas- 
ing for each iteration. We currently use h = 0.5 for the temperature steps and 
/j = 0.8 for the level population steps. 

The partial derivatives needed for Eq. l4.68l are computed numerically, with 
a step-size ofAx/x = 0.01. T here are metho ds for computing the optimal value 
of this step-size, discussed in lHeathl (119970 . 

The criteria for convergence are difficult to specify. We mark convergence 
when the relative Newton-Raphson step, |Ax,/;c,|, is less than 1% in all direc- 
tions. For the thermal equilibrium, we require AT /T < 1%, and also that the 
net heating rate is smaller than 1 % of the total heating rate. For the ionization 



'which is the most efficient technique for large systems. 
^Today's CPU's can do ~ 10 flops per clock cycle. 
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balance, we similarly require that the net ionization rate is smaller than 1 % of 
the total ionization rate. 

The excitation structure is more difficult to quantify a convergence criterion 
for, because there are often r elativel y large spurious changes for high levels 
with low populations. iLucy (l2003bl) uses a criteria that the relative change 
averaged over the five lowest levels should be smaller than 10^^. The exact 
choice of convergence criteria is still under development for the code here. 

Handling of run-away situations has also been implemented. Temperature 
steps leading to negative temperatures, or temperatures above 3 x lO'* K, are 
halved until they give a value within the — 3 x 10^ range. Similar boundary 
checks on the ionization and excitation solutions have been implemented. 
When convergence fails for some zone within 100 iterations, the zone is 
flagged, and the solutions are set to their previous ones. Final models are then 
checked for flags, and if any are present the error has to be tracked down by 
more detailed investigation. 
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The 20th century witnessed the development of radiative transfer theory and 
techniques for computing the radiation field inside and outside astrophysical 
objects. The basic theory is rather straight-forward, but the resulting equation 
systems are often large, complex, and difficult to solve. The fundamental dif- 
ficulty is that the radiation emitted at one point in the object may influence 
the physical conditions at another point, which emission may influence the 
conditions at the first point, etc. Much of the Uterature and research is focused 
on how one may construct approximate equations that are easier to solve than 
the exact ones, and how one solves the original equation system in the most 
efficient way. 



5.1 Basic theory 

The flow of radiative energy at frequency V, position x, time t, and in direction 
n is described by the specific intensity 

= dtdZvda ' ^"'^--^Hz-Ster- , (5.1) 

which is a complete description of the radiation field, apart from polarization. 
The mean intensity is the angle-averaged value of /y 

Jv{x,t,n) = (^Iy{x,t,n)dQ. . (5.2) 

The change in /y as the beam travels from point x to point x+dx comprises 

emission, represented by the emission coefficient jV (erg s^' cm~^ ster~^), 
and absorption, represented by the absorption coefficient tty (cm~^): 

Iv{x+dx,t +dt ,n) — Iy{x,t ,n) = jv{x,t ,n)ds — av{x,t ,n)Iy{x,t ,n)ds , (5.3) 

where we let ds be a distance along the direction n. 

The ratio of emission and absorption coefficients is called the source fiinc- 
tion 

Sy(x,t,n) = ^^^^ , erg s-^cm-^Hz-iSter-i , (5.4) 

av{x,t,n) 

which has the same dimension as /y. 
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Rewriting the LHS of Eq. l5.3l as 



Iv{x + dx,t + dt,n) —Iv(x,t,n) = 

= -:-^dt + -^ds 
dt ds 

dly dt dly 
= —ds + -z—ds 
dt ds OS 

Id d\ 
-cJt^JS)'^'''^ 



(5.5) 



we obtain 



]_d_ d_ 
c dt ds 



/v(x,f,n) = 7v(x,f,n) - av(x,f,n)/v(x,f,n) 



(5.6) 



In flows with V/c <^ 1 and no radioactivity, the time-derivative 
term is much smaller than the space derivative and may be ignored 
(IMihalas & Weibel Mihalaslll984l) . In SNe, the requirement for dropping the 
time-derivative is 

— < mm 1, ■-] . (5.7) 

c t 

In spherically symmetric models, we can replace x with radius r, and direc- 
tion n with the angle /i relative to the radial direction, so that 

/v(x,f,n) ^/v(r,?,/l) 
7v(x,f,n) jv{r,t,IJ.) 
av{x,t,n) — > av{r,t,}x) . 

In spherical coordinates, Eq. 15.11 (with the time-derivative dropped) be- 
comes 

:i-M^) d 



dr 



dpi 



Iy{r,t,n)=jy{r,t,n)-av{r,t,n)Iy{r,t,n) . (5.8) 



As we shall see, jv{r,t,iJ.) generally depends on /i -integrals of /y, and this 
is therefore a partial integro-differential equation in r and /i. Solutions may 
be approached by difference equations of the form in Eq. 15.81 or by integral 
equations. For moving media like SNe, the in tegral equation s become cum- 
bersome, and the differential form is preferred (IMihalasll 19781) . 



5.1.1 Emission coefficient 

The emission coefficient, jy{r,t,iJ.), is different for different observers due to 
Doppler effects (Sect. 15.21 ). Below, we derive expressions for the coefficient in 
the rest frame of the emitting gas. 

Emission of photons may occur by spontaneous processes, stimulated pro- 
cesses, or by scattering. Because stimulated rates are proportional to /y, just as 
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absorptions, they are most easily treated as negative absorptions (Sect. [5X21 ). 
so only spontaneous processes and scattering contribute to the emission coef- 
ficient. 

Spontaneous emission involves a change of energy state for the electron, 
which at any given time is either free (^kinetic + ^^potentiai > 0) or bound 
(£^kinetic + ^^potentiai < 0). Depending on the state of the electron before and 
after the photon emission, emission processes can then be divided into 
free-free (ff), free-bound (fb), and bound-bound (bb). In total we have 

jv=/j+ji'+fj'+rr' ■ (5.9) 

Free-free emission j{,^ 

A free electron can emit photons if it is accelerated in the Coulomb field of 
an ion. This is called free-free emission, or Bremsstrahlung. The emission 
coefficient for ion k (representing some element and ionization stage), is 

j^^^' = ^^-^^' ^^'^'^emZlg'jfT'/^ , (5.10) 

where is the electron number density, is the number density of the ion, 
and gjy^ is the Gaunt factor. Free-free emission will play a role mainly at early 
times when ionization and temperature are high. Just as for H II regions, it is 
not very competitive for radiating thermal energy. 

The total free-free emission coefficient is the sum over all ions 

7f = IjP. (5.11) 

k 

Free-bound emission /J' 

When a free electron is captured into a bound state / of an ion k by radiative re- 
combination, free-bound (recombination) emission occurs. The emission co- 
efficient for reactions with ion k+\is 

= -laf-*^'(r)«,«,+i/iv, v>Xki 

= 0, v<Xki, (5.12) 

where n^ is the electron number density, n^-^i is the ion number density, 
av^'*'(r) the reaction rate, and Xkl is the ionization potential. Often, the 
detailed v-dependence of the reactions rates al^'^' are unknown. In such 
instances, an approximate treatment is to distribute the total reaction rate 
a,{jf'^'(r) (which often is known) over a box-shaped profile of width 
Mv = ksT: 

jmi ^ _L^f^^„^„,^^/jv, xki <hv< (xki+hAv) 

= , otherwise . (5.13) 
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This approach has been used in this thesis for all atoms except hydrogen, 
where exact free-bound reaction rates are available. 

The total free-bound emissivity is obtained by adding up the contribution 
by all ions and levels 

;f =£;f'^'. (5.14) 

k,l 

Bound-bound emission j*^ 

Excited, bound electrons transition spontaneously to lower states by certain 
transition rates A„/ (s^'), where u and / denote the upper and lower states. The 
emission coefficient is 

where i//^(v — Vg"') is the emission line profile. In SNe, the velocity gradients 
are so large that the typical Doppler shifts become much larger than the ther- 
mal line widths. We can therefore approximate the line profile with the Dirac 
delta function. 

Va(v-v^"') = 5(v-Vo^"0 . (5.16) 

As shown in Sect. 15.2.21 in SNe one may replace the emissivity above with an 
effective emission coefficient 

.hbMi ^ J-n,Akuilikui{nki,ntu,t)hv5{v - v^"') , (5.17) 

471 

where Pkui is the escape probability of the line, and then skip doing the radia- 
tive transfer through the line itself. 

Bound-bound deexcitation can occur also by two-photon (2/) emission, 
where the sum of the photon energies equals the transition energy. The line 
profile is then a smooth function, starting at Vq"' and extending to infinity. In 
Type II SNe, two such transitions are important; H I 2^— ls{X = 1215 A — oo) 
and He I 2s — ls{X = 601 A — oo). The H I two-photon emission is stronger 
than Lya as long as 

(5.18) 



t < 8 years -—^ r .(5.19) 

\10MqJ V3000kms-V 

However, boosts in the escape probability due to overlapping lines may pro- 
duce other estimates. 

The total bound-bound emissivity is 

jtb ^ £ jthMi (5_20) 

k,uj 
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Scattering j^y"" 

Scattering refers to photon-matter interactions where the photon continues on 
(in a new direction) with only a small frequency change. In the code, scattering 
is not included in the computation of the emissivity function, but is instead 
treated directly in the radiative transfer simulation (Paper III, Appendix A). 
However, for completeness I describe the process for constructing scattering 
emissivity below. 

Scatterings in SNe mainly occur on free electrons, called Thomson scatter- 
ing. The emission coefficient is 

;T'(m) = ^Orn, <f^£ Iyi^')R{v',v,^' ,^)da'dv' , (5.21) 

where the cross section is Gt = 6.65 x 10^^^ cm^ (independent of wavelength), 
and R{v' , v , jJ.' , jJ.) is the redistribution function. If we ignore the small energy 
transfer occurring and take v = v', we can write 

;r"(M) = 4^6.65 X \0-^\ej^IvWg{^',^)dQ! , (5.22) 

where the angular redistribution function is 

g{n',H) = ^[l+cos\n'-H)] . (5.23) 

Another scattering process is Rayleigh scattering, which however has a much 
smaller cross section than Oj, and is neglegible in SN environments (see Sect. 
|5X2]l. 

Isotropy 

All spontaneous emission processes are isotropic in the rest frame of the gas 
since they involve thermal distributions of particles. Scattering, on the other 
hand, is anisotropic since it depends on the angular distribution of the incom- 
ing radiation field, and on the angular redistribution function How- 
ever, several_au^ors_^ndtha^his anisotropy introduces quite small effects 
(e.g. lAuer & van Blerkoml[l972l) . so one can to good approximation take 



= constant . (5.24) 
The total emission coefficient is then isotropic: 

jv{r,t,lJ.) jy{r,t) . (5.25) 



5.1.2 Absorption coefficient 

Just as for the emission coefficient, the absorption coefficient av{r,t,iJ.) is 
different for different observers. Below, we derive the expression for the coef- 
ficient in the rest frame of the absorbing gas. 
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Absorption may occur by free-free processes, bound-free processes (pho- 
toionization), bound-bound processes, or by scattering: 

ay = (4^ + a^f + a^b^ ^cat (5^6) 

ff 

Free-free absorption ay 

The free-free absorption coefficient is 

= 3.7 X 10«r-i/Vv-3 (1 - e-'^'l^'^^) Y,4nkg)f , (5.27) 

k 

where summation is done over all ions k. The free-free optical depth at 1 jJLm 
is (ignoring the exponential, which is small) 

^inm ""-'^ '5000 Ky^QM^j 1^3000kms-V \0.\ J VlOOdays 

(5:28) 

It may therefore influence the early nebular phase, in particular the infrared 
spectra. Because of free-free absorption, we can not expect to see any far- 
infrared emission lines early on. 

Bound-bound absorption 

The line absorption coefficient is, dropping the k index, 

«^^''" = <T/«<^>(v - v^)n/ - (7„/V^(v - v^)n„ , (5.29) 

where <^(v — Vg") is the absorption Une profile, and ^'^(v — Vq") is the emission 
line profile (which is the same for spontaneous and stimulated emission). The 
second term is the correction for stimulated emission. The cross sections are 
related to the Einstein B-coefficients as 

aiu = -^Bi^ , (5.30) 

and 

Oui = -^Bui . (5.31) 
An 

Using the Einstein relations Biu/Bui = gu/gi waAAui/Bui = IhvlJc^, and ap- 
proximating the line profile as a box function of width Av, we get 

«v = -^—^ulKloini n«)-r- , |V - V/„o| < -TT ■ (5.32) 

81 8u Z 

Since in homologous expansion dV jds = \/t, and the Doppler formula gives 
^lu = cAv/V/„o = '^vA/mO; a line profile of velocity width AV/„ corresponds 
to a spatial distance of 

Liu = mul^=^vXiuot , (5.33) 
as 
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hbju hh.lu 



which is called the Sobolev length. Then, assuming a constant value of oCy*''" 
over the distance L/„, we obtain the line optical depth 

Liu = ^'-AiAltm(l-^'^) , |v-Av|<^, (5.34) 
871 gi \gunij 2 

which is called the Sobolev optical depth (z^). The result holds for all line 
profiles. Note that since the density evolves as p{t) oc t^^ oc the Sobolev 
optical depth scales as (assuming constant ionization/excitation), just as 
in static media. 

The size of the line interaction region relative to the size of the nebula is 

Z^nebula Vt V V \AmJ ' ' ' ' 

At r = 5000 K, AViu is of order 10 km s^^ for hydrogen. For V = 10"* km s^^, 
we then get Lune/^nebuia ~ 10^^. It is thus a good approximation to assume 
constant physical conditions over the line interaction region as long as (1) 
clumping does not occur on a smaller scale than ~ 10^^, and (2) the line is not 
so optically thick that the Lorentz wings are important (then AVh, > /^yf'''""""'). 

h f 

Bound-free absorption 

Whereas calculation of bound-free absorpti on coefficients are a vailable for the 
ground-state of all atoms (we use data from Verner et al. 1996^. they are more 



scarce for excited states. For hydrogenic atoms, the coefficient scales as 

= 7.9 X gnZ^ , (5.36) 

where n is the principal quantum number. For non-hydrogenic ions, this ex- 
pression is still a good approximation for levels with effective quantum num- 
bers n > 10. For other levels, the cross sections are more complex, often with 
resonances at certain frequencies that may differ by several orders of mag- 
nitude from the underlying level. However, the effective 'smearing' of fre- 
quencies in the differentially expanding SN environment likely allows lower- 
resolution approximations to be used to good accuracy. We use Eq. 15.361 for 
all excited levels with unknown bound-free cross sections. 

Just as stimulated bound-bound emission is most easily treated as nega- 
tive absorption, stimulated recombinations are most easily treated as nega- 
tive bound-free absorptions. Then, the absorption coefficient from level / is 



( Mihalaslll970 ) 



V ni/n* J 
where n^ is the number density of the next ion, is the number density of 
the ionizing level, and starred quantities represent LTE values. At the relevant 
wavelengths and temperatures (~ 10^ A and ~ lO'* K), the exponential has 
value ~ 10^^. We therefore ignore stimulated recombinations in the transfer. 
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Scattering ay"" 

The scattering coefficient is, assuming contribution only from electron scat- 
tering 

^^catt ^ ^scatt ^ g^g ^ 10-25„^ . (5 38) 

The optical depth is, for a homogeneous hydrogen sphere 

= 2 8 ( ^ ' ('-^] ( ' 

■ \10MqJ V3000kms-i J \^o.01 J V 100 days 

It is indeed electron scattering that traps the radiation for the first ~ 100 days 
in Type II SNe, causing the photospheric phase. 

Another scattering process is Rayleigh scattering, where photons scatter 
on bound-bound transitions of much higher energy than the photon energy 
(Vui ^ v). The cross section is 

ar"'^ = ar/»/ ^ (5.40) 

where /„/ is the oscillator strength. The Rayleigh cross section is much smaller 
than Ot since v <C v„/ . However, if ^ 1 , the total absorption coefficient may 
come to exceed the Thomson one, as Lyman lines or metal resonance lines on 
neutral elements provide scattering targets. In low-velocity, low ^^Ni-mass 
SNe, Rayleigh scattering may not be neglegible in early nebular phases. 

Isotropy 

All absorption coefficients are isotropic in the rest-frame of the gas: 

av{r,t,n) ^ av{r,t) . (5.41) 

Fig. 15. II shows the various continuum optical depths at 150 days for a 10 
M0 hydrogen sphere with V = 3000 kms"\ T = 6000 K, = 0.05, and an 
LTE excitation structure. We see that Thomson scattering (dashed line) and 
photoionization (solid line) are still important for UVOIR photons, whereas 
also free-free absorption (dash-dotted line) is important in the FIR. 

Fig. 15.21 illustrates the line opacity by plotting the probability for a photon 
emitted in the core of the 19 M0 model in paper III to escape the SN, at 
300 days. We see that below ~ 5000 A there is complete blocking, between 
5000 — 12, 000 A there is partial blocking, and after that there is little blocking. 

5.1.3 Lambda iteration 

Clearly, the emission and absorption coefficients, jv and ay, depend on the 
state of the matter (temperature, ionization, excitation), and on ly, both di- 
rectly (scattering terms) and indirectly (through the influence of ly on the 
state of the matter). Even if LTE is assumed, this coupling remains because 
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Figure 5.1: The continuum optical depths for free-free absorption (dash-dotted), 
bound-free absorption (soHd), and Thomson scattering (dashed), for a 10 M0 hy- 
drogen sphere at 150 days, with V = 3000 kms"', T = 6000 K and Xe = 0.05. The 
bound-free absorption has been computed up to principal number n = 4, above which 
LTE becomes a bad approximation. 



the temperature is influenced by the radiation field. Thus, equation |5. 8 l is more 
difficult to solve than it may appear, and has to be solved in conjunction with 
the equations governing the the state of the matter (Chapter 3). 

If scattering is treated as isotropic, we in general have, for a 
time-independent problem, that jy and tty depend on 7v' at the same location, 
for all frequencies v'.The mean intensities in turn, depend on the emission 
and absorption coefficients at all other points and frequencies: 

Jy{r) = A[jy,{r'),aAr% Vv',Vr' , (5.42) 

where the K-operator represents computation of Jy (r) from jyi (/) and tty' ('"') 
by some method. 

To solve the coupled equation systems, one starts with a guess for iy{r) and 
ttv ('')■ A radiation field can then be computed by solving the radiative transfer 
equation (15.8b by some method. With this solution, new evaluations of the gas 
state can be made through the formulas described in Chapter 3. This allows 
new jy and ay to be computed, and the procedure is repeated. 

In the stellar atmospheres context, A-iteration has proven to work well in 
the case of small or moderate optical depths, but converges slowly or not at 
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all for large optical depths ( Mihalas 1970l) . The reason for the latter is that in- 
formation is propagated only one optical depth per iteration through the grid, 
so will, in a random walk, take ~ steps from boundary-to-boundary. In the 
thick case, one may still devise tricks to enhance the convergence properties, 
generally referred to as accelerated lambda iteration (ALI) or operator split- 
ting procedures (e. g Cannon 1973). The first widely used methods of this kind 



were developed by lScharmeii(ll98ll : ll984) . 



In the nebular-phase SN context, the A-iteration scheme has quite benign 
properties since the optical depths are moderate, and in addition the gas state is 
mainly determined by radioactivity rather than by the UVOIR radiation field. 
As we saw i Chapter 3, the radioactive ionization, excitation, and heating rates 
are only weakly coupled to the gas state. The coupling between conditions at 
different points in the nebula is thus greatly reduced, and the convergence 
properties are improved. It has been demonstrated, with the code in this the- 
sis and a similar code developed at MPA in Munich, that A-iteration without 
any strong ac celeration schemes is rapidly convergent in nebular phase SN 
environments (IMaurer et al.ll201 ih . 

Still, in the early phases radioactivity reaches only the inner parts of the 
nebula, and conditions in the outer parts are completely determined by the 
radiation field. I f optic al depths become high, ALI methods are needed. 



Eastman & Pintol ( 1993) developed an ALI operator sp littin g technique for 
this sit uation, and'Lucvl tm^ l2003bL Eooi) . Isiml (l2007h . andlKromer & Sim| 
(l2009h developed Monte Carlo methods enforcing radiative equilibrium by 
the use of indestructible energy packets. The radiative equilibrium condition 
is otherwise just fulfilled asymptotically, and the enforcement from the first 
iterations speeds up the convergence. 



5.2 Moving media 
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5.2 Moving media 

If the medium is moving, radiation transport is complicated by the fact that 
the absorption and emission coefficients differ between the rest frame of the 
gas (which is the one we have worked in so far) and the rest frame of the 
star (the observer frame), becoming direction-dependent due to the direction- 
dependence of Doppler shifts, aberration, and advection. Whereas aberration 
and advection involve effects of order V jc, Doppler shift s introdtice ma jor 
changes already when V/Vthermai becomes non-significant ('Mihalas 1978 ). It 
is therefore reasonable to take Doppler effects into account but neglect aber- 
ration and advection. 

The relationship between the observer-frame (unprimed) and comoving 
frame (primed) emission and absorption coefficien ts are (note that the pri mes 
were not written out in Sect. 15. 1.T1I5. 1.21 ) (e.g. Rybicki & Lightmanlll979h 



Jv(At) 



7. 



V 1-^^ COS /J 



V(r) 



cos \i 



(5.43) 



and 



av(At) = 



a 



V 1 



V(,-) 



COS \i 



y(r) 



COS \i 



(5.44) 



where the last equalities are valid for spherically symmetric flows only. The ef- 
fect of the motion on the continuum transfer is small, but for the line-transfer 
the situation may change completely if the frequency shifts become larger 
than the line widths. The angular dependencies makes solution of the radia- 
tive transfer equation more complicated, and dire ct solu tion in the observer's 
frame is only feasible in certain cases ( Mihalaslll978h . It is instead advan- 
tageous to compute all radiation-matter interaction in the comoving frame, 
the frame where the interacting matter is at rest, because here emission and 
scattering coefficients are isotropic. Since specific intensity can be shown to 
transform as 

U=l'y:(-^ . (5.45) 



V' 



One may thus solve the differential equation of transfer in the comoving frame 
to update /(,, and then use Eq. 15.451 to transform to the the lab frame or to 
the next comoving frame of interest. This is the method preferred in most 
computer codes, also used in this thesis, even though we transfer Monte Carlo 
packets rather than /y explicitly. (Sect. 153] ). 

Since W jc even in SN explosions, it is not a bad approximation to con- 
sider only the wavelength shifts and set = and = a^, ( Lucy 1999bl; 
20051). This is also justified by the moderate accuracy to which many other 
physical processes can be modelled in SNe, and the fact that we ignore the 
time-derivative term in the equation of transfer. 
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5.2.1 Scattering 

A fundamental difference to static media is that scattering now introduces po- 
tentially large frequency changes for the photon in the star's frame. Consider a 
photon that moves in the same direction as the scattering particle, and is scat- 
tered 180° backwards. Given that the scattering is coherent in the rest frame 
of the scattering particle, the photon suffers a redshift of Av/v = 2y/c in the 
star's frame, where V is the velocity of the scatterer. If the photon instead 
moves anti-parallel to the scatterer, there is an energy gain of 2V /c. Since the 
nebula is radiating, photons must find themselves moving outwards (same di- 
rection as the gas) more often than inwards (opposite direction as the gas). The 
net effect of scattering is therefore that energy is transferred from the radiation 
field to the gas. 

Electron scattering will tend to broaden and redshift the line profiles be- 
cause any photon that escapes will have gone through additional frequency- 
var ying direction changes, with more redshifting than blueshifting ones (see 



e.g. lAuer & van BlerkomI 119721) . The corresponding increase in path length 
for the photons increases the chance of line or continuum absorption, and so 
the escape of line radiation from the core will be made more difficult. The 
rest-frame broadening due to thermal motions can be seen for the narrow lines 
in some Type Iln SNe, but is generally unimportant for the ejecta lines, where 
the broadening is dominated by the bulk flow. 



5.2.2 Line radiation 

In static media, a photon of a given frequency faces absorption by one or 
a few lines close to that particular frequency. Absorption can happen at any 
place in the nebula. In a moving medium, on the other hand, each Une is only 
in resonance with the photon at one particular point along the trajectory. This 
means that a velocity gradient provides an escape mechanism for the photon, 
shifting most of the atoms out of resonance with it. 

On the other hand, there may now be other lines that shift into resonance 
with the photon somewhere along the path. Thus, each photon can now be 
absorbed by many lines, but for each line just at one point. The effect of many 
lines being coupled together is called interlocking, and their collective opacity 
is referred to as line blocking. 

Sobolev approximation 



SobolevI (119571) derived the solution for the line transfer in a moving medium 
with expansion velocities much larger than the intrinsic line widths. With the 
assumption of an isotropic and frequency-independent s ource functio n Sui, the 
mean intensity weighted over the line profile becomes (iCastorlll97ol) 

/„, = + ^7 = ( 1 - Sui + P!,Jt, , (5.46) 
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where ^'/"^ is the contribution to the mean intensity due to emission in the line 



itself, 7fjf is the contribution due to external photons, which approximately 
jb 
■'ill- 
escape probability, which for homologous flows is 



equals 7^^,, the mean intensity in the blue wing of the line, and j8^^, is the Sobolev 



ti 



1 



(5.47) 



ul 



where T^^^ is defined by Eq. 15.341 The escape probability jS^^ is the angle and 
frequency-averaged probability that a photon emitted in the line will escape 
rather than be reabsorbed in the same line by another atom. The net radiative 
rate between the levels is 



Mui = UuAui + UuBuiJui - niBuJui 
= UuAui + riuBui 



{l-P^,)Sui + l5^t4i 



■ mBiu 



nAulPul + i^uBui - niBi^)4il5^ 



where we have used 
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niBiu - fiuB, 



{l-^',i)Sui+P!A 

(5.48) 
(5.49) 



(5.50) 
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This means that instead of emitting a total flux «„Aj,/ and doing the radiative 
transfer through the line, we can emit a flux riuAui^f^i and start the transfer 
after the line. We also recover the result stated in Sect. l4.6.Tl that the photoex- 
citation rate is given by 

Rl^=BiM. (5.51) 
Generalized escape probability 

The actual source function does not look like Eq. 15.501 but rather like 
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(5.52) 



Hummer & Rybickil (Il985h and lChugail (Il987h derived expressions for the es- 
cape probability when the a^^^ term is non-zero, having a contribution from 
continuum opacity. They find that one can write, denoting R = Oi"^^ ja}^^'^ 



■"ul 



Pf,+R-f{Tli,R), 



(5.53) 



where f{x^^,R) is a function of order unity. For complete redistribution over 
the Doppler core 
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. (5.54) 
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For partial redistribution, IChugail (119871) finds 

'a\i/4 



/(t;^„/?) = 1.87 



R 



where a is the damping parameter. For T = 5000 K and Z = 1 : 



a = 5^<^=4xl0-4^^'<"^"' 



4kAv 



10^ s-^ I V 5000 A 



(5.55) 



(5.56) 




Figure 5.3: The boost in escape probability due to continuum destruction as function 
ofR^ a™"7a/j^ The solid curves are for complete redistribution and the dashed 
ones for partial redistribution, using a= 10^^. 

As an example, for Lya in the core hydrogen zone in our model for SN 
2004et at 400 days (Paper HI), R = 10"^^ /p, ~ 300, and ^'ff = 3 x 10"^° ~ 
4j8'^. In Fig. 15.31 we show the boost in escape probability as function of 7? for 
two different values of T"^. 

The escape probability may be further increased by absorption in neighbor- 
ing lines. This has particular importance for the Lya line. The n = 2 popula- 
tion in H I will be affected if AjSjf^ > Aiy, i.e. ^^l^ > 4 x 10"*^. Several metal 
lines exist around Lya that may bring ^21'^ above this value. Table l5.1] lists all 
lines within 2 A (w SOAV"*^"™' at 7 = 5000 K) on the red side of Lya having 
'^Lya ^ ^^^^ hydrogen zone in the 12 M0 progenitor model of 
Paper IV, at 300 days. It is clear from the Tfot/ zf^^ values that jSjf ^ may tran- 
scend 4 X 10^^. Preliminary simulations have given jS^;^^ 10^^ (Jerkstrand 
& Fransson, in prep.), which is the value used in paper HI and IV. 



5.3 Monte-Carlo method 
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Line 






X 


1 L\OL 




Fe II 


1215 


85 


A 


4.5 


2.7 X 10" 


-9 


Fell 


1216 


24 


A 


14 


1.2 X 10" 


-9 


Fell 


1216 


27 


A 


15 


1.2 X 10" 


-9 


Fell 


1216 


52 


A 


21 


1.3 X 10" 


-9 


Fell 


1216 


85 


A 


30 


1.2 X 10" 


-9 


CrII 


1217 


14 


A 


37 


8.3 X 10- 


10 


Fell 


1217 


15 


A 


37 


1.5 X 10" 


-9 


Fell 


1217 


62 


A 


49 


1.0 X 10" 


-9 



Table 5.1: Lines between 1215.67 and 1217.67 A, in the core hydrogen zone in the 12 
Mq model in Paper IV, at 300 days. The x values are the line separations in number 
of (hydrogen) Doppler widths at 5000 K. 



Another possibly important line overlap occurs between Lyj3 (1025.72 A) 
and O I 2p3p-3d3D (1025.76 A). The Lyj3 line-width is AA = 0.04 A at 5000 
K, so the line separation is within the Doppler core. One may show that the 
absorption probability in the O I line is 



oBowen 
PLy/3 



2.8 X 10 



-5 



/o 



fo. 



solar 



(5.57) 



where fo is the number fraction of O I relative to H I, and /o,soiar is the solar 
v alue. Thi s is the value we use in paper III and IV. 



Chugail (119801 demonstrated that the Sobolev approximation holds also 



under partial redistribution ( /' |)"^ / t^^n becomes frequency-dependent), 



which iHummer & Rybickil (11992 ') verified by numerical simulations. How- 
ever, they made the important point that the spatial random walk performed 
by the photon during the scattering process is likely to take it so far away 
from its starting poin t that one can hard ly as sume the phys ical conditions to 
stay constant. Earlier, lOsterbrockl (Il962h and Adarnsj (Il972h made the impor- 
tant findings that in static nebulae with optical depths T > 10^, photons escape 
by random- walking to the edge of the nebula rather than by diffusing to a fre- 
quency in the wing where the nebula is optically thin. This situation will occur 
also in SNe at some optical depth, and the Sobolev approximation will then 
break down. 



5.3 Monte-Carlo method 

The transport of photons (or any kind of particles) may be computed by fol- 
lowing them on their trajectories through the grid. The interactions along the 
way are determined by random draws, and the method is therefore referred 
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to as the Monte Carlo (MC) method. The technique is especially attractive in 
situations of multiple, complex interactions that are difficult to construct and 
solve the corresponding equation systems for. 

For photon transport, the emission and absorption coefficients jv and a^, 
derived in Sect. l5.1.T] and [5X2l can be used to construct probabilities for the 
interactions. There is no need to compute /y explicitly, but instead the emer- 
gent spectrum is constructed by simply binning the escaping photons in ac- 
cordance with their wavelengths, angles, and number (if photon packets are 
used). In Paper II we give a detailed account of our implementation of the 
creation and transport of photons using a MC technique. In summary, we cre- 
ate photons by random sampling of the emissivity function jv{r), send them 
off in a random direction, and let them be absorbed when the integrated op- 
tical depth (computed from the a{r) function) transcends a randomly drawn 
optical depth. By constantly transforming the photon frequency to the local 
comoving frame, we can use the rest-frame expressions for the absorption co- 
efficients. Following an absorption, random numbers determine a deexcitation 
channel, and new directions for the scattered/fluorescent photons. 

Apart from its conceptual simplicity, the MC method has the further 
advantages of being relatively easy to construct and test, suitable for 2D 
and 3D problems, easy to implement parallel algorithms for, and suitable 
for anisotropic phase functions as easily as isotropic ones. The numerical 
difficulties often encountered in attempting to calculate /y are to a large ex tent 
removed when photons are instead moved in MC schemes (lAueiil2003h . In 
the context of moving media, the MC method is especially suitable since 
the line interlocking (Sect. 15.2.21 ) produces long and complex expressions 
when attempting to solve the equation of transfer by integration. MC transfer, 
on the other hand, simply moves the photons from line to line in a scheme 
simple to code. 

Another advantage, exploited in this thesis, is the possibility to include 
clumping and mixing by direct incorporation into the random sampling ma- 
chinery for the photon trajectories. This is explained fully in paper II (Sect. 
2.1.1). 

The disadvantages of the MC method are that computing times can become 
long, and that rare channels are difficult to study due to statistical noise. If 
the optical depths are too high, the method becomes inefficient as the pho- 
tons have to be tracked through a large number of interactions. But with ever- 
growing computer speeds and increasing availability of multi-processor com- 
p uting, MC methods are likely to keep increasing in popularity. 



Abbott & Lucyl (Il985h demonstrated the use of the MC method for deal- 



ing with the expanding atmospheres of hot stars. The explosion of SN 1987 A 
prompted the development of MC codes for SNe (ILucvI 19871 iMazzali & LucM 



2005 



1993 ) . Some of these MC codes use indestructible energy packets (e.g. lLucy 



Kasen et al.l 120061 ISimI 120071) . so that after any radiation-matter inter- 



action, a packet with the same energy as the incoming packet is re-emitted. 
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This ensures that the radiation field satisfies the condition of radiative equi- 
librium (e.g. it is divergence-free), even if the temperature structure does not. 
This, in turn, has been shown to accelerate the convergence properties of the 
A-iteration (,Lucv.. 1999b: ,2002) . 



5.3.1 Probability distribution sampling 

The core of the MC method is to sample probability distributions. To do this, 
we need the cumulative probability distribution 

ia P{x')dx' 

where p{x) is the (not necessarily normalized) probability distribution over the 
variable interval [a,b\. By sampling from a uniform distribution z = [0, 1] 
and solving for x, we sample the p{x) distribution {the. fundamental princi- 
ple, e.g. Duderstadt and Martin 1979). For example, isotropic emission corre- 
sponds to p{0) = sin0, where 6 = [0,7r]. We get 

"^^^^-/o^inew 

1 — COS 6 

~ 2 

->cos0 = l-2z (5.59) 



Another example is the distance a photon travels before being absorbed. Here 
p{t) = e^"^, so 

^T = -ln(l-z) (5.61) 

Some probability distributions are not analytically integrable. We then ac- 
cumulate the integral numerically until the drawn z value is reached, 
s, but overestimated absolute fluxes. 



5.3.2 Line interaction 

It is the coupling between radiation field and line transitions that c auses the 



most difficulties in obtaining convergence for the A-iteration process (iMihalas 



It is therefore of interest to consider approximations in the radiation- 



line interaction, and here the MC method allows us great freedom in tailoring 
specific treatments. 

When a photon is absorbed in a line, several things may follow. The atom 
may de-excite in the same transition, re-emitting an identical photon in a new 
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(random) direction. Since there is no net energy transfer between matter and 
radiatio n, this type o f inte raction gives ben i gn con vergence properties. It was 



used by iLucyi (119871) and iMazzali & Lucyl (Il993h for computing the transfer 



tiirough SN atmospheres. 

The atom may also deexcite to another level with the emission of a fluo- 
rescent photon. This photon will be of higher energy than the original one if 
the fluorescence occurs to a state below the original level, and of smaller en- 
ergy otherwise. Here, demanding a divergence-free radiation field means that 



Lucv 


( 1999bl). 


Mazzali ( 


2OOOI) and 


Kasen et al. 



(l2006h implemented treatments allowing for the initial fluorescence channels 
by weighting the transition probabilities with the transition energies. This 
gives the correct relative fluxes in the fluorescence lines, but incorrect absolute 
fluxes. 

The atom may also transition to any other level by a collision with an elec- 
tron. The probability that a thermalizing deexcitation occurs is 



Pthen 



Lf~>thermal 
l<u^ul 



(5.62) 



h/J 



The collision coefficients c^'j'""""' are of order lO^^Ug s^^ If any allowed, 
optically thin channels exist, Auip^^j is of order 10^ s^\ and thermalization 
cannot occur since it would require «g > 10'"* cm^^, much higher than in 
SNe. If only forbidden channels exist, A„/j8^; is of order 10^^ — 1 s^\ and 
thermalization occurs if rie > 10^ — 10^ cm~^. This electron density is achieved 
in most parts of the ejecta for quite a long time. Consequently, radiative energy 
absorbed to meta-stable levels will thermalize. 

A lower degree of matter-radiation coupling can be achieved if the 
photoexcitation/de-excitation rates due to non-local radiation are ignored. 
One may t h en compute the excit a tion structur e by LTE approxi mation 
dLucvl Il98i iMazzali & Lucvl Il993[ ILucvI ll999bL iKasen et aimOO eS or by 
approximate NLTE solutions (paper II). Then, we must continue following 
the excited atom in the current iteration, since the absorption will not be 
remembered in the next one. For this thesis, such a decoupled scheme 
was used in paper I and II. Photon absorptions were then followed by 
randomly sampling radiative and coUisional deexcitations, until the atom had 
de-excited to the ground state. For papers III and IV, an upgrade to include 
photoexcitations for a limited number of transitions was implemented. 



5.4 Implementation 
5.4.1 Assumptions 

The major assumptions for the radiative transfer treatment are 



5.4 Implementation 
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Spherical symmetry. Although the choice of the MC technique is partially 
motivated by the planned use of multidimensional models in the future, 
only one-dimensional ones have so far been considered. There is growing 
evidence that strong asymmetries are present in SN ejecta, but for the Type 
IIP SNe analysed here, much of the action occurs in the large hydrogen 
envelope where asymmetries should play a smaller role. 
Sobolev approximation with no line overlaps (in the comoving frame) ex- 
cept for Lya and Lyj5. A detailed discussion about the possi ble problems 



associ ated with the Sobolev approximation can be found in iKasen et al. 



(l2006l) . Other uncertainties likely introduce errors more serious than the 
ones caused by this approximation, a statement supported by test calcula- 
tions (Eas tman & Pinto 1993) . Line overlap can be a problem in the photo- 
spheric phases when the ejecta are hot and dense, causing a huge number 
of lines to be optically thick. In the nebular phase, at lower temperature and 
density, it is less of a problem. 

No time delays (the photons travel instantaneously between points). See 
Sect. [33] for a discussion about this point. 

No effects of order V/c except for the frequency Doppler shifts. Relativistic 
effects may easily be included, but has not been implemented in considera- 
tion of neglecting time-delays and time-derivatives. For the SNe considered 
in this thesis, V/c is small enough for these effects to be neglegible. 



5.4.2 Transport of free-bound continua 

A common method in astrophysics is to replace the recombination rate a by 
a — ags, and assume that all free-bound continua to the ground state are reab- 
sorbed on-the-spot. The justification is that these continua usually have very 
high optical depths. Here, we perform transfer also for these continua, as they 
may affect the photoionization rates of trace elements as well as be partially 
absorbed by lines. 



5.4.3 Packet splitting 

In runs where photons are tracked over their whole trajectories, many MC 
codes avoid splitting the photon upon matter interaction, but instead choose 
one of the resultant ph otons by random draws. This is a matter of coding 
simplicity (lLucyll2002h . However, we find that photon- splitting is not a big 
problem, and has the advantage that fluorescence is sampled to a higher degree 
of accuracy. At high enough optical depths, though, the memory requirements 
for following chains of fluorescence photons become demanding. 
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5.4.4 The randomized grid 

A major idea advanced in this thesis (paper II) is to exploit the nature of the 
MC transfer by letting clumping and mixing be represented by a randomized 
ejecta structure. Few objects are more suitable for such a treatment than SNe, 
where the ejecta consist of several chemically distinct clumps that become 
mixed around during the explosion. To represent such a structure in a fixed 
one-dimensional grid, each zone type inevitable has to be represented by a set 
of shells, completely enclosing the material inside. This may give emergent 
spectra that differ drastically from spectra from clumpy models. One example 
is the impossibiUty to have some optically thick clumps while still allowing 
for some clear line-of-sights through the ejecta. 

5.4.5 Spatial resolution 

We assume constant conditions in each zone type within the core. The uniform 
density and ^^Ni distribution implies that the radiation field is not expected to 
change by more than a factor ~2 over the core. 

For the envelope, we wish to resolve changes by a factor of ~2 in density or 
radiation field. The density changes by a power law no steeper than p (r) ~ 
in most parts of the hydrogen envelope. For a r^" power law, the required 
resolution is 

— =2'/", (5.63) 
ri 

which equals 1.2 for n = 4. This is the value we use. 



5.4.6 Line-list 



To reduce run times, it is desirable to consider only lines with optical depth 
above some threshold value in the MC simulation. This value should be signif- 
icantly smaller than 1 , since many weak lines can produce collective blocking. 
We typically use T^j^^ = 10~^. The ty pical number of lines in the early nebular 



phase is then 10 — 10 . iLucyl (Il987h uses a cut-off at r 



.5 

mill 



10 



-2 



5.4.7 Wavelength resolution and range 

To optimize run-times, as narrow wavelength range as possible should be used 
for the radiation transport. We certainly want to include the He I ionization 
edge at 504 A. Little emissivity occurs at shorter wavelengths than this, so a 
lower limit of Amin = 400 A is reasonable. The choice of the upper limit has 
to be at least the wavelength where continuum absorption may occur. Since 
the Balmer continuum is optically thick for up to a few years, = 3650 
A is a reasonable choice. However, at early times both Paschen and free-free 
absorption may occur, and A^ax should then be at least 8000 A. 
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With lower limits Amin and Amax> and wavelength resolution AA /A = r, the 
number of frequency points is 

log(l+r) r V^in/ 

It is thus more rewarding to set r as large as possible than to squeeze the 
frequency range. We certainly want to have r <^V /do resolve the lines well. 
Since V/c ~ 1 x 10"^, r ~ 1 x lO"^ 

is a good choice. 

In the UV, lines are separated by less than this in the early phases. However, 
the requirements of the Sobolev approximation are anyway not fulfilled here, 
so retaining this approximation while using higher resolution in the transfer is 
not meaningful. 



5.4.8 Code testing 

An extensive library of testing routines was built up alongside the develop- 
ment of the code. This library has verified the basic functionality for all units, 
and allows an opportunity to redo the testing after any code changes. 

In addition to this, the output of the code was compared to the output of a 
similar code simult aneously developed at MPA. The results are described in 



Maureret all ( 2011 ). 
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6 Summary of the papers 



6.1 Paper I: The 3-D structure of SN 1987 A's inner 
ejecta 

The explosion mechanism of supernovae is still unknown (Sect. I3.1I ). The 
main aim of this paper is to diagnose the explosion of SN 1987A by determin- 
ing the three-dimensional structure of the ejecta. To this end, SN 1987A was 
observed in the near-infrared with the VLT SINFONI integral-field spectrom- 
eter, which provides a spectrum for each spatial location in the field of view of 
the instrument. SINFONI is an adaptive optics assisted instrument and there- 
fore operates in the near-infrared. Identifying strong lines in the spectrum, and 
assuming that there is no blending with continuum or other lines, the velocity 
of the gas can then be determined, giving the full 3D structure of the ejecta. 

The results revealed an asymmetric morphology, elongated more towards 
the equatorial ring plane than perpendicular to it. This probably rules out any 
rotationally-driven explosion mechanism, as one then would expect the latter 
type of asymmetry. The morphology instead favors the SASI-type of explo- 
sion, described in Sect. 13.1.21 because this type of instability causes asymme- 
tries in random directions. One may note that a random direction will, with 
50% probability, be within ±30% of the equatorial plane, producing an elon- 
gation apparently quite parallel to the equatorial plane, as observed here. 

My most important contribution to this project was to run spectral modeling 
to determine the origin of the 1.644 /xm and 2.058 /xm lines, i.e. to identify 
which elements and which zones emit them. Both lines turned out to have 
a surprising origin. According to the model, the 1.644 /im line is not domi- 
nated by the [Fe II] 1.644 /im transition, but by [Si I] 1.645 /xm. This line is 
emitted by silicon in the Si/S zone, and tells us that the 3D mapping shows 
the distribution of oxygen-burning ashes in the SN (see Sect. 12. 1 .5] and l3. iTb . 
The theoretical predictions for the amount of silicon produced (~ 0.1 M0), 
and its mixing with radioactive nuclides (with ~10% of the ^^Ti), are in good 
agreement with the observations. 

The iron distribution is instead traced by the He 1 2.058 ^um mapping, as the 
model shows that this line is produced by helium created by freeze-out in the 
iron core, described in Sect. 13.1.11 This conclusion is further strengthened by 
the fact that the He 1 2.058 /xm and [Si I] 1.645 jxra distributions are similar, an 
unlikely outcome if the He I emission comes from the helium zone and not the 
Fe/He zone. The results therefore provide confirmation for a-rich freeze-out 
actually having occured in the explosion. 
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For the rest of the spectrum, a reasonable model agreement was obtained 
using MC^'^Ti) = 1 x 10~^ and and optical depth for the dust in the ejecta 
"Z^dust = 1> although below ~ 1.4 /im, the model is too low by a factor of ~ 2. 
This can likely be explained by the higher ^^Ti-mass found in Paper II (1 .5 x 
10^^ M0), and possibly also by input by X-ray s at late times (Sect . I4.2.31) . 
The X-rays do not penetrate to the core material / Larsson et al. 2011 ), which 
means that they have no direct influence on the [Si I] 1.646 /im and He 1 2.058 
Hm lines, can can thus not explain the observed morphology. However, the X- 
rays may influence the H, He, and continuum emission from the envelope, or 
fluorescence thereof. 



6.2 Paper II: The ^^Ti-powered spectrum of SN 1987A 

This paper contains a lot of the basic background for the modeling, and results 
of application to the late-time optical spectrum of SN 1987A. The spectrum 
of SN 1987A in the ^^Co and ^'Co-powered phases have been thoroughly 
modelled before. Less attention has been paid to the '^^Ti-dominated phase 
(f > 5 years). At these late times, the ejecta temperatures have dropped to 
a few hundred degrees, and the entire UV/optical/NIR spectrum is produced 
by non-thermal processes. The powering is mainly due to positrons, as most 
of the 7-rays escape the ejecta (see Sect. 14.21 ). These epochs therefore offer 
a unique opportunity to see emission from the deepest interiors where the 
positrons get trapped. 

We choose to analyze the epoch of eight years, where a high-quality HST 
spectrum as well as a ground-based NIR spectrum exist. Previous attempts to 
identify lines in the HST spectrum has le ft question-mark s for several of the 
strongest lines, especially in the UV (Ch ugai et al.lll997 l). With the addition 
of multi-line radiative transfer to our model, we were in a unique position to 
identify these lines. After eight years, increasing contamination by the cir- 
cumstellar interaction makes spectral analysis of the ejecta more difficult. 

An important question is where in the core the positrons deposit their en- 
ergy. As described in Sect. 14.21 even a weak magnetic field will trap the 
positrons close to the radioactive material. If no magnetic field is present, 
they can instead travel substantial distances, but will still be absorbed some- 
where in the core. We find a model without a magnetic field to produce too 
strong emission lines from the oxygen and hydrogen components, and there- 
fore conclude that a magnetic field is indeed present. Energy deposition into 
the mainly neutral iron clumps produces prominent Fe I emission lines in the 
spectrum, matching several of the previously unidentified lines. 

From investigation of the line formation process, we could show that multi- 
line transfer plays an important role, even at eight years after explosion. As 
described in Sect. l5.1.2l the velocity gradients in SNe lead to very efficient line 
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blocking, and a large fraction of the emergent spectrum is made of scattered 
light. 

The dust clumps that formed at ~ 500 days (Sect. 13.6.21 ) are still optically 
thick at eight years, damping the red sides of the emission lines. With the 
best estimate for the effect of the dust on the UVOIR field, we determined the 
^"^Ti-mass to I -S^q j x 10^^ M0. Interestingly, this is close to the value in Cas 
A (1.6lo| ^ M©)' the only SN where the y-ray decay lines by '^'^Ti have 
been directly detected. As we discuss in the paper, the determined "^^Ti mass is 
higher than obtained in spherically symmetric models, but enhanced produc- 
tion occurs in asymmetric explosions. The results here therefore ties in with 
the findings in Paper I, that the SN 1987A explosion was quite asymmetric. 



6.3 Paper III: The progenitor mass of the Type IIP 
supernova 2004et from late-time spectral modeling 

In Chapter 2, we describe how the amount of nucleosynthesis occurring in a 
star is strongly dependent on its main-sequence mass. The aim of this paper is 
to constrain the amount of nucleosynthesis products present in the bright Type 
IIP SN 2004et, thereby providing constraints on its progenitor mass. This SN 
is interesting because the progenitor was detected in pre-explosion images, 
and analysis of these suggests a significantly less massive star (14 M©) than 
obtained by hydrodynamical models of the SN light-curve (> 25 Mq). With 
an excellent UV to mid-infrared nebular phase data-set, spectral modeling can 
add a crucial third diagnosis for which type of star SN 2004et was. 

To analyze the nebular-phase spectra, we apply the spectral synthesis 
code, finding the nebular phase spectrum to best agree with models of a 
~15 M0 (non-rotating) progenitor (as described in Sect. I2.4[ rotation would 
result in the progenitor mass being lower for a given amount of metals 
produced, so 15 Mq can be considered an upper limit). We identify Mg I] 
4571 A, Na I 5890, 5896 A, [O I] 6300, 6364 A, blends of metal Unes at 
1.15 — 1.2 /im and 1.4-1.6 /im, and [Ne II] 12.81 /^m, as the most useful 
lines for constraining nucleosynthesis in Type IIP SNe. We can therefore 
conclude that our spectral modeling is in good agreement with the progenitor 
analysis, but in disagreement with the hydrodynamical modeling. 

We can also confirm that silicate dust emission, as well as SiO and CO fun- 



damen tal band emission, is present in the Spitzer spectra taken by lKotak et al, 



(l2009h . SN 1987A did not have any significant dust continuum until after ~ 
500 days, whereas in SN 2004et it is present already at day 350. The CO emis- 
sion is in agreement with the total heating of the CO zone in the 15 Mq model. 
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6.4 Paper IV: Constraining the physical properties of 
Type IIP supernovae using nebular phase spectra 

This paper looks at a sample of Type IIP SNe with progenitor detections as 
well as spectroscopic data sets in the nebular phase, investigating the variety 
in line profiles and flux ratios. How large is the diversity within this class of 
SNe, and what can the differences tell us? Standard stellar evolution models 
predict all stars in the 8 — 30 M0 range to end as red supergiants and explode 
as Type IIP SNe. But all the progenitors dete cted in the last deca de or two tend 
to cluster in a smaller range of 8 — 17 M0 (ISmartt et al.ll2009h . Can spectral 
analysis settle which of the ranges is the correct one? 

A first step is to understand the formation of the strongest lines in the neb- 
ular phase spectra. As most SNe lack infrared spectra, we concentrated on the 
optical [O I] 6300, 6364 A , Ha, [Fe II] 7155, 7172 A, [Ca II] 7291, 7323 A, 
and Ca 11 IR triplet lines. As most of the SNe in the sample have identified 
progenitors between 8 — 15 M0, we interpret the evolution of various lines by 
the use of a 12 Mq model. 

For the ca l cium lines, we could show that the conclusion drawn by 



Li & McCrayi (119931) for SN 1987 A, that they originate in primordial calcium 



in the hydrogen envelope, is valid for all epochs. From the discussion in 
Sect. 12. 2[ it is clear that the final hydrogen envelope mass varies little with 
the progenitor mass, and the envelope lines therefore show little variation 
within the sample. Furthermore, none of the SNe in the sample show any 
flat-topped emission lines, which means that the hydrogen-zone material has 
been mixed down towards zero velocity. This occurs by the hydrodynamical 
instabilities described in Sect. 13.1.21 and the results therefore suggest that 
these instabilities invariably occur when the star has a significant hydrogen 
envelope left at the time of explosion. 

We also found that the [O I] A A 6300, 6364 doublet has more than half 
of its flux coming from natal oxygen in the hydrogen envelope, a possibility 
not previously recognized. This fact helps to explain why the [O I] AA6300, 
6364 to Ha line ratio does not vary much within the sample, even though the 
synthesized oxygen mass varies greatly with progenitor mass. 

My most important contributions to the paper consists of modeling the line 
flux ratios (Sect. 5), as well as modeling SN 2008bk to identify the emission 
lines (Sect. 6), but I was also involved in the development and formulation of 
all other chapters. 
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Publications not included in this thesis 



I NERO - A Post Maximum Supernova Radiation Transport 

Code 

Maurer I., Jerkstrand A., Mazzali, P.A., Taubenberger S., 
Hachinger S., Kromer M., Sim S., Accepted for publication in 
MAWA5.ArXiv: 1105.3049. 



n 



X-ray illumination of the ejecta of SN 1987A 

Larsson J., [3 authors], Jerkstrand A, [23 authors], 201 1, Nature, 
AlA, 484 
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